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Abstract 


In this thesis, we study the d —orthogonal polynomials and the d— orthogonality in the 
Sobolev sens. We présent an algebraic theory of classical d—orthogonal polynomials and 
we want to hll in some gaps. We broaden and close sonie inclusions that exist and are 
known perhaps as conséquences. Several characterizations of classical d-OPS are given in 
ternis of d + 2 terni récurrence relation as wcll as in ternis of functional and differential 
équations. An additional tool in determining the intégral représentation for such class is 
presented. 

Furthermore, we study the inner products which generate the sequence of orthogonal 
polynomials in the sense of Sobolev. 

Using the theory of Riordan group and d— orthogonal polynomials, we shall show that 
sequences of d— orthogonal polynomials, can be also generated by Riordan group. We 
interpret sonie families of d— orthogonal polynomials in the frame of Riordan group. 

keyword: 

Orthogonal polynomials , d—Orthogonal polynomials , classical polynomials , Riordan 
Matrix 
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Résumé 


Le thème principal de ce travail porte sur les polynômes d— orthogonaux au sens de 
Sobolev. Après avoir cité quelques propriétés de l’orthogonalité de Sobolev, nous avons 
construit des produits Scalaires qui génèrent des suites d— orthogonaux au sens de Sobolev 
et nous avons montré que les suites 2—classiques par rapport au vecteur r de dimension 2— 
des formes linéaires ne sont pas definie une suite 2— classique au sens de la généralisation 
de Sobolev 

Nous avons ensuite étudié de nombreuses caractérisations de la d —SPO classique. 
Tout d’abord, nous donnons une nouvelle démonstration de l’équation de Pearson en util¬ 
isant l’orthogonalité faible ainsi que sur une nouvelle caractérisation du d—SPO classique 
en ternies de la d— quasi-orthogonalité. D’autres caractérisations ont également été ex¬ 
primées en termes de la d—quasi orthogonalité ainsi que de l’orthogonalité faible. En 
outre, nous présentons un outil pour déterminer la représentation intégrale de telle classe 
de polynômes. 

Le second thème dans cette thèse est consacré au polynômes d—orthogonaux et les 
groupes de Riordan. Nous avons pu montrer que les suites de polynômes d—orthogonaux 
peuvent être également générées par ce type de groupe. De plus, nous avons exhibé 
quelques familles d—orthogonaux à partir des coefficients de ces matrices. Le résultat 
principal obtenu dans cette étude est que l’inverse d’une matrice de Riordan exponen¬ 
tielle L est un tableau des coefficients d’une famille de polynômes d—orthogonaux si et 
seulement si la matrice Stieltjes est égale à la matrice Production. 

Mots clé: polynômes orthogonaux , d—Orthogonalité , polynômes classique , Matrices 
de Riordan 
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Chapter 

Introduction 


1.1 Orthogonality 

The notion of orthogonal System of fonctions appeared through the study of certain 
problems of fonctional analysis (intégral équations, strum-Liouville problcm and more 
generally, boundary problems in partial differential équations). Orthogonal polynomi- 
als in general, and classical orthogonal polynomials in spécifie, hâve been the subject of 
extensive work. They are related to many problems of applied mathematics, theoretical 
physics, chcmistry, theory of approximation, and several other mathematical branches. In 
particular, their applications are widely used in théories such as Padé approximants, con¬ 
tinuons fractions, spectral study of Schrôdinger, discrète operators, polynomial solutions 
of second-order differential équations, and others. 

The first family of orthogonal polynomials known as ” Legendre polynomials” appeared 
with the first work of A.M.Legendre on the planetary movements in 1784 . Actually, lie 
established sonie common properties and lie derived a second order differential équation 
which lias as solutions this families. He also studied its zéros. 

Subsequently, other families were introduced. For instance, Hcrmite (1864), con- 
structed a new family, called Hcrmite polynomials, which was used in the interpolation 
theory. 

N. H. Abel and V. L. Lagrange P .L. Chebyshev began with T.J Stieltjes the création 
and development of general theory of orthogonal polynomials by the use of continuons 
fractions. 

E. N .Laguerre introduced in 1897 a family of polynomials (called Laguerre polynomi- 
als) and showed the link between it and the continuons fractions. 

In the middle of the nineteenth century, Jacobi introduced a new family which general- 
ized the polynomials of Legendre and Chebyshev. These three families (Jacobi, Laguerre 
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and Hermite), are known as the ”classical polynomials”. 

Other authors such as T. J. Stieltjes, L. M. Humbureger, T. Carleman, F. Hausdorff, 
and M .H. Stone hâve shown that orthogonal polynomials solve also the problem of mo¬ 
ments via continnous fractions. 

Bochner in 1929, demonstrated that the stnrm -Liouville differential équation, 
dehned by 

a(x)y" + fi(x)y’ + 7 (x)y(x) = A y(x) (1.1) 


has a polynomial solution if and only if dega(x) < 2 , deg/3(x) < 1 and y( x ) = 0 . 
Furthermore, by using a linear change of variable, Bochner has also shown that there 
are only 4 families of classical polynomials orthogonal solutions of (1.1) which are Jacobi, 
Bessel, Laguerre and Hermite. 


The second characterization of classical orthogonal polynomials is the Rodrigues 
formula which gives us the explicit expression of the classical orthogonal polynomials in 


ternis of a and f3 coefficients in équation (1.1). 


In 1938, Krall demonstrated that the classical polynomials satisûed in fact an even 
order differential équation of the type 


2 r 


a,(x)y (,) = -':</■ 


( 1 . 2 ) 


i= 1 


The results of Bochner and Krall hâve been generalized by Louriero and Maroni in 


2008 and 2010 by giving the link between the coefficients of ( 1 . 1 ) and ( 1 . 2 ). 


Other criteria for the classification of the classical orthogonal polynomials hâve been 
recently given by Kil. H. Kwon and L. L. Littlejohn (1996) on a Sobolev class, i.e., a 
classification based on a symmetric bilinear form dehned by 


q) = (ffi P q) + (b pV) , p,q£V, 

In 2006, Jamei and Koepf gave a class of a differential équation which generalizes the 
Bochner équation and the Sturm-Liouville problem. More precisely, they considered the 
case when the degree of the coefficients a and (3 are > 2 and also when they are rational 
functions. They showed that the solutions of this class are symmetric and orthogonal and 
they gave explicit expressions of their weight function. 

The notion of orthogonality has undergone a great deal of generalization. It began 
with the notion of (1 /p), (p > 1) orthogonality (A.Boukhmis 1988), the d-orthogonality 
(P.Maroni 1989), the biotogonality (Brezinski 1992), and ûnally the multiple orthogonal¬ 
ity (Aptekarev, van Assche ...). These generalizations intended to make the phenomena 
modeled by differential équations more compréhensible, to improve the approximation 
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methods, and to increase the order of the differential équations having these polynomials 
as solution. 

The notion of d —orthogonality appeared in the thesis of J.Van Isegham [82] . Tliis 
notion was used in the study of the Padé approximation of simultaneous formai d—sériés. 
Subsequently, it was observed that d—orthogonal polynomials satisûed an récurrence re¬ 
lation of order (d + 1). 

After two years, the first properties of a sequence of d— orthogonal polynomials hâve 
been given in 0. The author gave a new algebraic approach showing many interesting 
characterization for the d—orthogonality, based on the orthogonality’s vectorial form. In 
addition, lie introduced the notion of d—quasi orthogonality in the same article. 

Many authors hâve attempted to improve tliis theory. Although the classical char- 
acter lias been approached by Douak and Maroni in numerous articles, other researchers 
(Y. Ben Cheikh, N. Ben Romdhane, K. Douak, I. Lamiri, A. Bokhemis, E. Zerouki, A. 
Saib, etc.) studied certain problems wliich lead to the construction of many d-analogues 
of the classical families and to the discovery of new ones. 

In 2006, A. Boukhemis and E. Zerouki found an interesting resuit where they deter- 
mined a linear differential équation which admits the 2-orthogonal classical polynomials 
as solution given by 

Rn,i{x)Pnlz{x) + Rn i3 (x)pj l 2 } :3 (x) + R n , 2 {x) P^%(x) + R n ^(x) P n+3 (x) = 0, 
with Rn ti (x) for i — 1, ...,4 are polynomials of degree less than or equal to i and 

RnA X ) = F n,l( x )S 3 (x), 

where S 3 (x) is a polynomial which does not dépend on n and whose degree is less than 3. 

The classical character, i.e. the sequence of orthogonal polynomials in which the 
derived sequence is orthogonal (O P S). Tliis is known as the property of Hahn. 

1.2 Riordan array 

The concept of Riordan arrays has been introduced in 1991 by Shapiro et al [73] with 
the aim of generalizing the concept of Renewal Array dehned by Rogers and hâve 
pointed out its connection to the Umbral calculus [68]. Their basic idea was to dehne a 
class of infinité lower triangular arrays with properties analogous to those of the Pascal’s 
triangle whose éléments are the binomial coefficients. 

The Riordan groups are particularly important in studying combinatorial identities and 
sums. Tliis concept has been investigated by Sprugnoli ca who pointed out the relevance 
of these matrices from the theoretical and practical point of view. In addition, his work 
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verified that many combinatorial sums can be solved by transforming the generating 
functions. So, we can see that the Riordan array concept is particnlarly important. For 
example, identities involving Stirling nnmbers can not be treated by methods related 
to hypergeometric functions. In |76], Sprugnoli paid attention to the identities of Abel 
and Gould. Further applications in reciprocal functions are discussed in [221 [83] , on 
subgroups of the Riordan group in Peart and Woan [66], on sonie characterizations of 
Riordan matrices in Rogers [(TF. and Merlini et al [üD] . and in many interesting related 
results in Cheon et al. ÏÏÜ, He et al [38], Nkwanta 021. Shapiro ra, Barry [7] and so 
fort h. 

Successively, some other aspects of the theory hâve been studied [ÜD]. Further gener- 
alizations of these last characterizations are given in [3ÜJ [3ü]J. 

The Shefïer sequence is a very general concept that includes many polynomial sequences 
as spécial cases. The concept of Shefïer Riordan group of Shefïer-type polynomials is 
dehned by He et al. [38]. Tliey prove the isomorphism between the Shefïer group and 
the Riordan group and présent an équivalence between the Riordan array pair and the 
generalized Stirling number pair. 

The lower triangular matrices and matrix factorization problems hâve catalyzed many 
investigations in recent years. Peart and Woodson [M] did some researches on this prob- 
lem. They showed that some classical Riordan arrays hâve triple faetorizations of the 
form L = PCF , where P, C, F are also Riordan arrays. We notice also that the inverse 
of the Stieltjes matrix represents a coefficient of the sequence polynomials HEIEHEL]. 
In [79], the authors présent new faetorizations, and provide extensive examples of families 
of classical polynomials. 

The main results of this thesis 

This thesis contains six chapters. 

Chapter 3 

In Chapter 3, we are interested in the study of orthogonal polynomials of Sobolev. To 
this aim, we first investigate properties of symmetric bilinear forms and in particular the 
one which is dehned by 


0(p, q) = (ffi p q) + (t, p'q ), p,q eV. 

We give necessary and sufficient conditions for orthogonal polynomials with respect to 
the symmetric bilinear form 0 given above, which satisfy the Bochner’s second - order 
difïerential équation. A classification of the orthogonal polynomial sequences of Sobolev 
solutions of the Bochner difïerential équation will be investigated also. 
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Finally, we study the regularity of symmetric d—dimensional vector forms 0 = 
(0o, 0i, -, 0d-i) T defined by 

0r(p,9) := (aS(x - c),p q) + ( T r ,p'q ') 

More precisely, we show in the case when, a = 0 and r is 2-classical functional 
moment, the sequence of classical 2-orthogonal polynomials satisûes a non-homogeneous 
third-differential équation of Boukhemis type. 

Chapter 4. 

Many of the properties of the classical orthogonal polynomials hâve nice extensions 
in this d-orthogonality setting: there will be a liigher order linear récurrence relation 
[TÜ1154] , there are nice functional équation (Pearson équation) (23U3D], linear combination 
(hnite-type relation) between the sequence and its dérivatives [58] and differential équation 
[19] , In addition, if the same ideas are to be found, each aspect of the theory is more 
complicated and the links between the different aspects are not so clear. 

This chapter deals with algebraic aspects of classical d-orthogonal polynomials. In 
the first section, we recall some débilitions and characterizations of the d —orthogonality 
which we need in the sequel. In addition, it is well known that there is exactly 2 d sets 
of d—symmetric d—orthogonal polynomials. We discuss liere again this resuit from the 
combinatorics point of view. The main results are listed in the tliird section. Many 
characterizations of the classical d-OPS are presented. First, we give a new proof of the 
Pearson équation in the algebraic aspects with the aid of the weak orthogonality as well 
as on a new characterization of the classical d-OPS in terms of d—quasi-orthogonality. 
Further characterizations also expressed in terms of d-quasi-orthogonality as well as on 
the weak orthogonality. Further, we présent a set of tools hclp to determinate the intégral 
représentation for a such class of polynomials. 

Chapter 5 

In this chapter, we shall show that the sequence of d-orthogonal polynomials, i.e., a 
sequence of polynomial satisfying a linear récurrence relation of d + 2- terms, can be 
also generated by a Riordan array. General référencés for orthogonal and d—orthogonal 
polynomials can be found in [241 54] , Links between Riordan arrays and orthogonal poly¬ 
nomials hâve been explored in OSIEZ!- Note that there are some examples in the literature 
of Riordan arrays that generate some families of polynomials sequence satisfying a récur¬ 
rence relation in four terms and more as well as a récurrence with variable coefficients 
[?, [6, [23] • In this chapter, two examples of d-orthogonal polynomials are presented, for 
which an ordinary Riordan array exists. The first one is the d-Chebyshev polynomials 
of second kind. The second one présents in fact the co-recursive d-Chebyshev of second 
kind. The next sections deal with Stieltjes matrix. We generalize the resuit of Peart and 
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Barry where we show that the corresponding Stieltjes matrix characterizes the sequence 
of d-orthogonal polynomials and we obtain the explicit form of that Riordan array. For 
the exponential Riordan array, we show again that Stieltjes matrix characterizes the se¬ 
quence of d-orthogonal polynomials. Some examples are discussed. Finally, we improve 
the Shefïer Riordan array by giving an example. 


1.2. Riordan array 



Chapter 



Preliminaries and notations 


The présent chapter reminds sonie concepts about orthogonal polynomials of a real vari¬ 
able and their fundamental properties. 

2.1 Orthogonality in general 

Let V be the linear space of complex polynomials of one variable and V its topological 
dual space. We dénoté by (u, P) the action of u and P G V . 

2.1.1 Free and monic polynomials 

Let {P n } n>0 be a sequence of polynomials with coefficients in C. 

Définition 2.1 The sequence {P n } n>0 is said free if and only if deg(P„) = n, Vn > 0 

Définition 2.2 We say that a free sequence {P n } n >Q is monic if each polynomial P n is 
written in the form : 

n—1 

P n (x) = x n + ^2 a n -kX n " k for n> 0 
fc =0 

Proposition 2.1 If {P n } n>0 is monic, then there exist a unique sequence {/3 n } and a 
unique array Xn,vi 0 < v < n, such that 

Po{x) = 1, Pi(x) = x - Po 

P n + 2 (x ) = (x - /3)p n + 1 - E”=0 Xn,vP v (x) , n > 0 


9 
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Définition 2.3 Let {P n } n>0 be a sequence of free polynomials, the function G(x] t) which 
can be developed into a power sériés of t is called the generating function of the sequence 
{P n } n>0 , if it is put in the form : 

G(x,t) = ^ c n Pn(x)t n , where c n 0 , Vrt > 0 . 

n> 0 


2.1.2 Regular and definite positive forms 

Définition 2.4 We call the dual sequence of the free sequence {P n } n>0 the sequence of 
linear forms {«„} defined by 

^n^Prn) ■ {^m Pm) &nmi TXjTXl ^ 0 , 


where is the duality hook. The linear form uq is called the canonical form of the 
sequence {P n } n > 0 - 

Proposition 2.2 i) J£$For each monic sequence, the dual sequence exists and is unique. 
ii) In this case, we hâve 


f3 n = (u n ,xP n (x)), n > 0, 

Xn,v = {u v , xP n+ 1 (x)) , 0 < V < n. 


( 2 . 1 ) 


Définition 2.5 The form £ is called regular (or admissible, or quasi-definite) if there 
exists a free sequence {P n } n>0 such that 

£(P n Pm) = 0,n^m, 

£{P 2 n ) ± 0 , n > 0 . 

Such a sequence is called orthogonal (or regularly orthogonal) relative to £. Such a 
sequence is free and unique up to a multiplicative constant. 


Lemma 2.1 For any £ E V' and any integer p > 1 ; the following statements are équiv¬ 
alent [54\ 1 55]/ : 

I) {£, Pp-i) 7 ^ 0, {£, P n ) = 0 ,n>p, 

II) 3À„ G C, 0 < v < p — 1 , A p _i 7 ^ 0 such that 


E P —1 

A v £ 


v =0 


Lemma 2.2 The sequence {P n } n>0 is orthogonal with respect to £. Necessarily, £ = 
\£o, A 7 ^ 0. In this case, we hâve Ï5R [55]/ 


£n={(£o,Pn)) 1 PnUo, U> 0. 


when £ is regular, let A be a polynomial such that A£ = 0. Then A — 0. 
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Remark 2.1 For any linear form £ and for any polynomial n, we can define the two 
forms D£ = £' and n£ by 


{£',/):=-{£, f), {*£,f):={£,irf), f e V. 
and thus, the dual sequence {ü n } n>0 °f {Qn}, is given by 

ü' n = ~(ji + 1) u n+ 1 , n > 0, (2.3) 

where Q n (x) = ( n + l)^ 1 P' n+1 (x) , n > 0. 

Définition 2.6 The sequence {£ n } n>0 defined by 

£ n := £(x n ), n > 0 

is called a sequence of moments of the form £. We say that the form /i is real if the 
moments £ n , n G N, are real. 


Theorem 2.1 

defined by 

are ail non null. 


49lThe form £ is regular if and only if the Hankel déterminants 
A n = det (£ i+j )? J=0 ,n > 0 


Proposition 2.3 [Efi The following conditions are ail équivalent 

1) {Pn} n > 0 orthogonal polynomials sequence relative to p. 

2 ) 


£(ir(x)P n (x)) 


0 if deg 7r < n, Vn > 1 
fi 0 if deg 7 t = n 


3) 


£(x m P n (x )) = A n ô nm , A n fi 0, o < m < n 


Définition 2.7 A form £ is said positive definite if (<£,7 t(x)) > 0 for any polynomial 
7T = 0 such that tt ( x ) > 0 for ail real x. 


Proposition 2.4 [24\ \Jfi^ If a form £ is positive definite, then it is regular and ail mo¬ 
ments of £ are real. 


Theorem 2.2 {2f\ (fiPj A form £ is positive definite if and only if it is real and the 
déterminants of Hankel A n are strictly positive for ail n > 0. 


2.1. Orthogonality in general 
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2.1.3 Fundamental récurrence. 

Theorem 2.3 [2fl \f9)Let £ be an admissible linear form on V and {P n } n>0 the orthogo¬ 
nal monic sequence relative to £. Then this sequence satisfies a linear récurrence of order 
two, i.e. there exist two sériés of numbers {« n } B>0 and {/3 n } n > 0 suc h that 


P 0 {x) = 1, Pi (x) = x - f3 0 

P n +2(x) = (x - f3 n+1 )P n+1 - a n+l P n (x) , n > 0 


with the regularity condition a n ^ 0 for ail n > 0 . 
b n x n ~ 1 + ... for ail n > 1 . we hâve, /3 n = b n+2 — b n+ 1 
(by convention). 


Moreover, if we set P n (x ) 


and a„ = 


^n+1 A n _i 

X 


where A 


= x n + 
-1 = 1 


Let {-Pn } n > 0 be a monic sequence of polynomials satisfying the recur- 
2.4 where a n 7 ^ 0 for any n > 0 .Then there exists a single linear form £ such that 


Theorem 2.4 

rence 


£(Po) = «0 and £(P m (x)P n (x)) = 0 for any n 7 ^ m. 


Corollary 2.1 Let {P n } n>0 be a monic sequence of polynomials satisfying the récurrence 


2.4. Then the form /i is positive definite if and only and a n > 0 for any n > 0. 


Définition 2.8 Let {P n } n >o be a monic sequence of polynomials satisfying the récur¬ 
rence 2.4 We call associated polynomial sequence of {P n } n> 0 the sequence denoted by 
{P n (x, c)} n>0 , for any c G M, defined by : 


P 0 (x, C ) = 1, Pi O, c) = X - Pc 

P n+2 0, c) = (x - f3 n+c+1 )P n+1 (x, c ) - a n+c+ iP n (x, c ) , n > 0 

Theorem 2.5 J ~24f/ Let {P n } n>0 be a sequence of orthogonal polynomials with respect to 
the positive definite form p, then each polynomial P n admits n real and simple roots for 
ail n > 1. Moreover, the zéros of P n and P n+ i are alternate, i.e. 


X 71 ■ 1 .v ^ Xn,v ^ ■ l 'n - 1 - 1 ^ ^ ^ l,îA 1, 

where {x n , v }v = i and {x n +i w}”=i are respectively the zéros of P n and P n+ \. 


2.1.4 Classical sequences 

Définition 2.9 Let {P n } n>0 be a sequence of orthogonal polynomials. This sequence is 
called classical if these sequences of successive dérivatives are a.lso orthogonal (Hahn’s 
property). 


2.1. Orthogonality in general 
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Theorem 2.6 


The classical sequences are the only real polynomial sequences 


solution of the second-order differential équation of the hypergeometric type of Sturm - 
Liouville 

a (x) y" + (3 (, x) y' + A n y = 0, (2.5) 


where a and (3 are polynomials of degree two and one respectively and where \ n is the 
following constant: 

X n = a"(x)n(n — 1) + n/3'(x ) 


Theorem 2.7 [2f}j The classical orthogonal polynomials are given by the formula of Ro¬ 
driguez 

P n (x) = 1 -^( a n (x)w(x)) : (2.6) 

k n w[x) dx n 

where k n is a non-zero constant and w is a solution of the differential équation ( aw)' = (3w. 


Corollary 2.2 If £ admits the following intégral représentation 


£(P) = / p(x)w(x)dx , for any polynomial P, 


then the sequences of polynomials given by 2.6 are orthogonal to the weight function w 


2.1.5 Integrated Functional Représentation 

We know, from Duran’s classical theorem concerning the problem of moments of Stieltjes, 
that every moment function can be represented by a measure in the form of an intégral 
of Rieman-Stieltjes. Now it is interesting to establish the intégral représentations of the 
moment functionnelle. 


Définition 2.10 [24] Let 0 be a function defined on / CI. 0 is said to be a rapidly 
decreasing function on / if, for every pair of integers we hâve 

„ d m 


x 




< Too for any x G I. 


dx m 

Définition 2.11 A bounded non -decreasing 0 whose moments 

/»+oO 


(2.7) 


Rn ■ 


x n d(f(x). 


are ail finite is called a distribution function, 

Theorem 2.8 fiTlff For any sequence {p n } n>0 of complex numbers there exists a fast decay 
function on R + (respectively on such that 


r»+OQ 




x n (p(x)dx respectively /i n := 


P+OO 


x n (f){x)dx. 


( 2 . 8 ) 


2.1. Orthogonality in general 
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Theorem 2.9 [fflj Let {-Pn} n>0 be a regular orthogonal polynomial sequence satisfying 
the differential équation 2.5. If w(x) is an orthogonalization weight function of {P n } n>0 
then w(x) satisfies the following differential équation 

(a(x)w(x))' — f3(x)w(x) = g(x) (2.9) 


where g(x) is a weight distribution with ail its zéro moments, i.e. 

(g(x),x n ) = 0,n > 0 . ( 2 . 10 ) 

Conversely, ifw(x ) is a distribution satisfying the following conditions 

I) w(x) decreases rapidly when |x| tends to oo; then ( w(x), x n ) exists and is finite for ail 
n > 0; 


II) w(x) defines a non trivial moment function, i.e., its moments are not ail null ; 


III) w(x) is a solution of équation \ 2.1 0\ in the distribution space in 
then w(x) is said orthogonalization distribution of {P n } n>0 . 


Définition 2.12 Relation 2.9 is called the non-homogeneous weight équation for the 


differential équation 2.5 


2.2 The d— Orthogonality 

Let us consider linear forms r 0 , ...,rd_i (d > 1). 


Définition 2.13 [53] The sequence {P n } n>0 is called d-orthogonal polynomials se¬ 
quence, in short a d-OPS, with respect to the d—dimensional vector of linear forms 

r = (r 0 ,...,r rf -i) T , if it fuimis 


(r a , x m P n (x)) = 0, n > md + a + 1, m > 0, 
(T'a, x m P md+a (x)) ^ 0 , m > 0 , 
for each 0 < a < d — 1. 

In this case, the d-dimensional functional — is called regular. 


( 2 . 11 ) 


Remark 2.2 [JS]1 If {P n } n is a d—orthogonal polynomial sequence, then its polynomials 


are exactly of degree n and can hence be normalized , thus it follows the uniqueness of 
the sequence. 

2) We say that the functional T is regular of dimension d 

3) When d = 1, we find the usual notion of regular orthogonality 


2.2. The d— Orthogonality 
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Theorem 2.10 [ôffllt is necessary and suffïcient, for the functional T of dimension d to 
be regular, that the following déterminants be non-zero : 


H , 


md-\-v 


1 

O 

Ti 

Y 

L md+v— 1 

Ti 

r 2 

r md+p 

T m—1 

T m 

^777,— \-\-TYld-\-V —1 

r 1 

m 

r 1 

1 m+1 

r 1 

(m+l)d-\-v— 1 

1 m 

pu 

1 m+l 

p v 

' ' ' (m+l)d-\-v—l 


, m > 0, v — 0, d — 1 


( 2 . 12 ) 


We hâve posed r n = r(a: n ) = (T 1 ^"), T 2 (x n ), ...,T d (x n )) T ,n > 0. Represents actually 
d line. 

Therefore, each déterminant H md+V is of order md + v if md + v > 1. 


By convention Hq — 1. Each of the ûrst m lines of the déterminant |2.12] represents 
in fact d line, 
we hâve (25) : 

m— 1 d v 

H md+V =n n u a (x u p ud+a - i(x)) n u^p^^x)) 

u= 0 a= 1 et —1 

for any m > 1 and 1 < v < d — 1 


m— 1 d 


H ™d =n n n m ^ 1 

n=0 a =1 ct=l 

v 

H v = \\U a {P a . l {x)) ,1 < v < d 


a=l 


and 


r“(x ri P n rf +Q _ 1 (a;)) = r f nd+n ,n >0,1 <a<d 

-H-nd+ot—l 

Let us now recall sorne characterizations we need in the sequel. 

Theorem 2.11 J/ Let {Pn} n >o be a monic sequence of polynomials, then the following 
statements are équivalent. 

(a) The sequence {P n } n>0 is d-OPS with respect to T = (r 0 ,..., Td_i) T , 


2.2. The d— Orthogonality 
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(2.14) 


(b) The sequence {P n } n>0 satisfies a récurrence relation of order d + 1 (d > 1 ): 

d-1 

Pm+d+l (x) = (x - (3 m+d ) P m +d (x) - ^2 'y^+d-vPm+d-l-u (x) , Tïl> 0 , (2.13) 

v=0 

with the initial data 

P 0 ( x ) = 1, Pi (x) = x- p 0 , 

P m ( X ) = (x - Pm-\) Pm-l (x) - J2P=0 'Jm-ï-uPm-2-u (x) , 2<Ul< d, 
and the regularity conditions 7^ +1 7 ^ 0 , m > 0 . 

(c) For each (n,is), n > 0, 0 < v < d — 1, there exist d polynomials 4>n U , 0 < /i < d — 1 
such that 

d—1 

«nd+v = Pn,u U ^ n > °) 0 < 1 / < d - 1, (2.15) 

fj,=0 

and verifying 

deg^ y = n, 0 < v < d — 1 , and if d > 2 , 

degç^ < n, 0 < n < v — 1 , ifl<z/<d — 1 , (2-16) 

deg^j, < n — 1 , î/ + l</i<d— 1 , if 0 < z/< d — 2 . 

Notice that the coefficients {/3 m }and { 7 m+1 }are given by 


p v = (u v , xP v ), 0 < v < d — 1 
7n+l+« = {Un+v, xP n+ d) , 0<V<d-l,7l>0 


(2.17) 


Proposition 2.5 j69i] For any d—OPS satisfying the récurrence relations 2.13 and 2.14 
the first moments 

(uo) î = Po 

(w 0 )i = Po + 7i 1 

(«1)2 = Po + Pi 

(^2)3 — Po + Pi + P2 

(« 3)4 = Po + Pi + P2 + Pz 

(« 1)3 = Po + Pl + ^2 + A)Al + +7l 1 + 72 1 

(“ 0)3 = Po + 7i 1 (2/?o + Pi) + 7i 2 

(« 2)4 = Po + Pl + P2 P P0P1 + P0P2 + P1P2 + 7i 1 + 72 1 + 73 1 

(«1)4 = Aj + ^1 + (Po + Pi)(PoPi + 27i x ) + 72 1 (/ 3 o + 2 /?i + / 3 2 ) + 7 i “ + 72 2 

(« 1)4 = Po + 7i 1 [ 2 / 5 q + (A) + Al) 2 + 7i 1 + 72 X ] + 7i 2 (2A) + A + A) + 7i ' 


2.2. The d —Orthogonality 
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Proposition 2.6 [56\ 1 The statements in theorem 2.11 are équivalent to the following 


(d) Xn+d- i, v = 0, n > v + 1 and Xv+d-i,v ^ 0, v > 0, 

(e) 

d— 1 

XU n U n —\ T “I - Xn+'U,Ti^'n+l+i') 71 ^ 0 (w_i 0) 

^=0 (2.18) 

Xn+d-l.n 7 ^ 0, n > 0. 

Corollary 2.3 fT^The multiplicity of zéros of any d — OP S is at most d. Moreover, 
any d + 1 consecutive polynomials as well as an d + 1 consecutive polynomials of 
r—associated sequence {P n } n>0 hâve no common zéro. And any r > 0 the polynomi¬ 
als Pn\ Pn +l \ .., Pn +l ^ hâve no common zéro 


2.3 Quasi - monomiality 

We give some results of quasi - monomiality in order to get the dual sequence of a given 
polynomial set by using the lowering operator associated with the involved polynomials. 
The results obtained will be applied to Boas -Buck polynomial sets. For an integer j, we 
dénoté by the space of operator acting on analytic function that augment (resp-reduce 
) the degree of every polynomial by exactly j if j G N (resp j G Z“), Z“ being the set of 
non positive integers.This tacitly includes the fact that, if cr G A^ 1 ), then cr(l) = 0. 

Définition 2.14 [g] Let a G A* h A polynomial set {P n } n>0 is called a sequence of 
basic polynomials for a if : 

(i) P 0 (x) = 1 , 

(ii) P n (0) = 0 whenever n > 0, 

(üi) crP n (x) = nP n -i(x). 

As a conséquence of this définition, we mention the orthogonality relation 

cr m P n (0) = n\S nm , n, m = 0,1,... (2.19) 

Starting from a polynomial set, We give a method of constructing its dual sequence in 
ternis of its lowering operator 

Lemma 2.3 Every a G A* b h as a sequence of basic polynomials 


2.3. Quasi - monomiality 
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Lemma 2.4 fEf Let {Q n } n >o b e a polynomial set and let a be its lowering operator. Let 
{P n } n>0 be a sequence of basic polynomials for a. Then there exists a power sériés 


such that 


ip{t) = y^ a n t n ,a 0 ^ 0, 

k=0 


<p{cr){Qn) = P n , n = 0 , 1 ,... 

Let us combine the two preceding lemmas, gives us the following theorem 


Theorem 2.12 JEfLet {Q n } n>0 be a polynomial set and let a lowering operator , and let 
{£ n } n> Q be its dual sequence. Then there exists a power sériés 


such that 


<p(t) = ^a n t n ,a 0 ^ 0, 

k=0 


(£n, f) = ^ )(/)(®)] x=0 = a ^ /(o), n = 0, 1, / G V 


( 2 . 20 ) 


If we apply Theorem 2.12 we obtain a further method to construct the sequence 
{£ n } n >o for general case from wliich we deduce the following expansion theorem 


Theorem 2.13 JE/ Let {Q n } n>0 be a polynomial set and let a be its lowering operator , 
and let {£ n } n>0 be its dual sequence. Then there exists a power sériés 


^ a n t n ,a 0 ^ 0 

k =0 

such that every analytic function f has the expansion 


/w = E 


0 ) 


n=0 


n\ 


Qn(^) 


( 2 . 21 ) 


If, morever, the translation operator T a commutâtes with a, then 


s ^ cr n tp(a)f(a) 

f(z + a) = 2^ -- Qn{z) 


n\ 


n =0 

Corollary 2.4 JEf Let {Q n } n> o be a Boas -Buck polynomial set generated by the formai 
relation 


G(x,t ) = A(t)B(xC(t )) = - t' 

n =0 ^ 


( 2 . 22 ) 


2.3. Quasi - monomiality 
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where 


A(t ) = B(t ) = ^ b k t k and C(t ) = ^ Cfct fc+1 (2.23) 

/c=0 /c=0 k= 0 

are i/ïree formai power sériés with the condition a 0 c 0 b k ^ 0, for ail k. Let u u x G 
A*” 1 ) such that 

oB(xt) = tB(xt) 

Put a = C*[y) where C*is the inverse of C i.e. 


C*(C(t )) = C(C*(t)) = t, with C*(t) = JV fc f fc+1 ,c* ^ 0, 

k =0 

then every analytic function has the expansion 

/w -t 

n=0 v ' 

If, moreover, the translation operator T a commutâtes with a, then 

'<‘ + ‘ 

n =0 v ' 

Définition 2.15 [9] A polynomial set {P n } n >o is called quasi -monomial if and only if it 
possible to define two operators a and r, independent of n, such 


crP n (x) = nP n -i(x) (2.24) 

rP n (x) = P n+1 (x ) 

Theorem 2.14 fVÿ Every polynomial set is quasi-monomial 

Theorem 2.15 Let (a, t) G A*- ^ x A^ 1 ) and let P := {Pn} n>0 be a polynomial set 
generated by 

a(x,t) = ( Z2S ) 

71=0 

Then we hâve the équivalences: 

crG(x, t ) = tG(x, t ) aP n = nP n _i 

and 
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where a := a x 

Définition 2.16 


tG(x, t) = ~G(x,t ) rP n = P n+l 

and r := r x 

[?] the operator a is called lowering operator of {P n (x)} n>0 


Remark 2.3 According to Theorems 2.14 and 2.15 , we deduce that if a function of 


two variables G(x,t) has an expansion of type 2.25, when degP n = n, there exist two 
operators (a, r) G A^b x Ab) such that G(x,t) is a eigenfunction of a := a^associated 
with the eigenvalue t and is also a solution of that équation associated with the operator 


2.4 Riordan arrays and orthogonal polynomials 

Riordan arrays give us an intuitive method of solving combinatorial problems, and liclp 
to build an understanding of many number patterns. They provide an effective method 
of proving combinatorial identifies and solving numerical puzzles as in [BTj rather than 
using computer based approaches [621 EU 

In this section we review mainly known results related to integer sequences and Riordan 
arrays that will be referred to in the rest of the work. 

2.4.1 Integer sequences and generating functions 

Définition 2.17 The ordinary generating function ( o.g.f .) of a sequence a n is the formai 
power sériés 

OO 

d( x ) = Ë a nX n (2.26) 

n =0 

Définition 2.18 The exponential generating function ( e.g.f .) of a sequence a n is the 
formai power sériés 

h( x ) = V'a n ^r (2.27) 

n\ 

n =0 

Définition 2.19 The bivariate generating functions ( b.g.f.fs ), either ordinary or expo¬ 
nential of an array (a n ,k) are d ie formai power sériés in two variables defined by 
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a(x,y) = y^ j a n ,kX n y k (o,g,f) and 


(2.28) 


n,k 



n,k 


Notation 2.1 The coefficient of x n is denoted by \x n \g(x) and from the définition of the 



l 


Lagrange inversion [35] provides a simple method to calcnlate the coefficients of the 


sériés reversion. 

Theorem 2.16 ÏÏÏBj Let <f{u) = J2™ =0 4>kU k be a power sériés of C[[w]] with 0 O fi 1 0. 
Then, the équation y = zfi(y) admits a unique solution in C[[n]] whose coefficients are 
given by 



n =0 


The Lagrange Inversion Theorem may be written as 



The simplest case is that of G(x) = x, in which we get 



If the product of two power sériés / and g is 1 then / and g are termed reciprocal 
sequences and satisfy the following. For o.g.fJs we hâve [161] 

Définition 2.20 A reciprocal sériés gfx) = Y^=o a n% n with a 0 = 1, of a sériés f(x) = 
b n x n with b 0 — 1, is a power sériés where g(x)f(x) = 1, which can be calculated as 
follows 


OO 


oo n 



(2.29) 


and for e.g.f.’s we hâve 
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^ ' ®rî I ^ ^ ^ ^ bi&n—i . i ®0 1 (2.30) 

ni ni 

n =0 n =0 

2.4.2 The Riordan group 

Définition 2.22 A Riordan array is a pair ( g(x ), f(x)) in which g(x ) and f(x) are formai 
power sériés such that g( 0 ) ^ 0 and /(O) = 0 . 

The pair defines an infinité, lower triangular array L = (l n ,k)n,km where 

L,k = [x n ]g(x)(f(x)) k 

If /'(0) 7 ^ 0 the Riordan array is called proper. 

Remark 2.4 From this définition, it easily follows that g(x)f(x) k is the generating func- 
tion of colnmn k in array 

Remark 2.5 A lower triangular matrix L is a Riordan Array if the generating function 
k th column is 

g( x )f(x) k , for k = 0,1,2,3,... 

with 


g(x) = 1 +g 3 x + g 2 x 2 + g 3 x 3 + 
f(x) = x + f 2 x 2 + f 3 x 3 + ... 

Theorem 2.17 The fundamental theorem of Riordan arrays 


«0 


bo 

CL\ 


h 








(9(z),f(z)) 


where the generating function of the two column vectors are A(z) = Y a n^ n and B(z ) = 

n> o 

Y b n z n respectively. The identity is true if and only if the following équation holds : 

n> 0 

9{z)A{f{z) = B(z ) 

Proof. We look at the Riordan Array (g(z), f(z)) column by column, and multiply it by 
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the column vector on the left-hand side 




ÜQ 



CL\ 

g gf gf 2 ■ ■ 


Û2 





This yields 


a 0 g + aigf + a 2 gf 2 + ... = g(a 0 + a x f + a 2 f 2 + ...) = g{z).A(f(z )) = B(z) 
and we hâve our resuit m 

Définition 2.23 The Riordan group 


(g(z), f(z)) / (g(z), f{z )) is a Riordan array and f(x) = f 0 + f x z + f 2 z 2 + ... , 
where / 0 = 0 , f x = 1 

i.e., 

each member of TZ is a lower triangular matrix with l's on the main diagonal.The 
multiplication in 1Z is 



(g(z),f(z))(G(z),F(z)) = (g(z)G(f(z)),F(f(z))). 
The identity is I = (1,^). 


The inverse of (g(z), f(z)) is 
of f(z),i.e. 


(- 1 


\9(f(z) 


, f(z) ), where f(z) is the compositional inverse 


/( f(z)) = f(M = z 
To check the inverse property, we compute 


g(f(z) 


, f(z) ) (g(z), f(z)) = 




g(f(z) ,/( f{z))) = ( 1 , 2 ) 


2.4.3 The A —sequence 

The A— sequence introduced by Rogers (1978) characterizes the column éléments of after 
the first column 
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Theorem 2.18 An infinité lower array D = (d nt k) n ,ke n 0 is a riordan array if and only if 
a sequence A = {a n } ngNo existe such that for every n, k G N 0 it is true 

n—k 

dn+l,k-\-l ^ ^ CLidn^k+i 
i=0 

even more, if D = (d(x), h(x)), then the generating function of the A—sequence is such 
that A{x) = a î> T? satisfies the équation 

h(x) = xA(h(x )) =>- A(x) = — 

h[x) 

2.4.4 The Z—sequence 

The Z—sequence introduced by Merlini et al characterizes the éléments of the first column 
of a proper Riordan as follows 

Theorem 2.19 Let D = (d(x),h(x)) = (d n ,fc)n,fce n- Then a unique sequence Z = 
((zo, zi, Z 2 , ■■■) can be determined such that every element in column 0 excluding the élé¬ 
ment in the first row can be expressed as a linear combination of ail the elementin the 
preceeding row with the coefficients identified as the éléments of the sequence Z satisfying 
the relation 


dn+ 1,0 ^ ^ Zidn,i G No 

i =0 

the generating function of the Z—sequence satisfies the équation 

<1{X) = 1 -x(Z(h(x)) ^ Z( ' r) = h(x) ( X “ d{h{x)) 

2.4.5 Stieltjes matrix 

In the context of Riordan arrays, We see the Stieltjes expansion Theorem in [65], defined 
as follows. 

Définition 2.24 Let L = ( l n ,k)n,k be a lower triangular matrix with l iti = 1 for ail i > 0. 
The Stieltjes matrix Sl associated with L is given by Sl = L~ l L where L is obtained 
from L by deleting the first row of L, that is, the element in the n th row and k th column 
of L is given by l Utk = l n+1}k 

Using the définition of the Stieltjes matrix above [63] leads to the following theorem 
relating the Riordan matrix to a Hankel matrix with a particular décomposition 
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Theorem 2.20 165 V Let H = = a n+ k} nk>0 be the Hankel matrix generated by the 

sequence 1, ai, a 2 , a 3 ,... . Assume that H = LDU where 


L (J”n,k) 


n,k>0 


1 0 0 0 0 


1,0 

1 

0 

0 

0 

'2,0 

h,i 

î 

0 

0 

-3,0 

h,i 

h ,2 

1 

0 

■4,0 

U,i 

li,2 

h,3 

1 


D = 


do 0 0 0 
0 d 1 0 0 
0 0 d 2 0 
0 0 0 d 3 


0, U = L t 


Then the Stieltjes matrix Sl is tridiagonal, with the form 


A) 

î 

0 

0 ■■■" 

«i 

A 

1 

0 

0 

OL 2 

A 

1 ••• 

0 

0 

« 3 

A» ••• 

0 

0 

0 

Oî4 






where 

A) = ai, «i = d\, f3k = h+i,k ~ h.k+i, ttfc+i = —r—, k > 0 

dk 

Now, we state other relevant results from this section, relating to generating functions 
which satisfy particular Stieltjes matrices. The first resuit relates to o.g.f.s 


Theorem 2.21 fhîj / Let H be the Hankel matrix generating by the sequence 1, ai, a 2 , a 3 ,... 
and let H = LDL T . Then Sl has form 


ai 

1 

0 

0 •••" 

«i 

A 

1 

0 ■■■ 

0 

a 

A 

1 ••• 

0 

0 

a 

A ••• 

0 

0 

0 

a • • • 
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if and only if the o.g.f. g{x) of the sequence 1 , 01 , 02 , 03 ,.. is given by 

9(x) = - --- 7 

1 — ai'X — aixj 

where 

f(x) = x(l + fif + af 2 ) 

Solving both équations above give us the required resuit. Similarly for e.g.f.s we hâve 
the following resuit. 

Theorem 2.22 [65fLet H be the Hankel matrix generating by the sequence 1, ai, 02 , 03 ,... 
and let H = LDL T . Then Sl has form 


A) 

1 

0 

0 ■■■' 

«i 

Ai 

1 

0 ■■■ 

0 

«2 

& 

1 ••• 

0 

0 

Oi 3 

As ••• 

0 

0 

0 

Q?4 * • * 






if and only if the o.g.f. g{x) of the sequence 1 , ai, a 2 , a 3 , ..is given by 

g(x) = J( a i — otif)dx, with g( 0 ) = 1 
where f is the solution of 

— 1 + Pf + a/ 2 , /(O) = 0 

The proof again in |65j involves looking at the form of the n th column of the Riordan 
array. However, intuitively tliis resuit can be from looking at the form of the matrix 
équation L = LS. In the case that L = [g(x), f(x)\ is an exponential Riordan array, we 
hâve the following 

Proposition 2.7 /?/ L = —L 

dx 

Proof. 


d 

dx 



OO 

n= 1 


j,n— 1 

(n — 1 )! 


'y ^ 9n +1 


71=0 


x ,L 

n\ 
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Equating the first column of matrix Land LS we hâve 

4~(g(x) = Po a( x ) + «i 9 (x)f(x) 
ax 

and second column equate to 

f ( x ) = Pif(x) + a 2 g(x)f(x) 2 

which gives the required resuit. ■ 

The Sticltjes matrix as we hâve seen above is a tridiagonal infinité matrix which is 
associated with orthogonal polynomials, and therefore hâve a generalization of the Stieltjes 
matrix to the Riordan group. Referred to as a production matrix ,it is deûned in the 
following ternis 

Let P be an infinité matrix (most often it will hâve integer entries). Letting ro be 
through throw vector r 0 = (1, 0,0, 0,we deûne r\ = r^iP, where i > 1. 

Stacking these leads to another infinité matrix which we dénoté by Ap .Tlien P is said 
production matrix for Ap . If we let U 1 = (1, 0, 0, 0,...) then we hâve 


A P — 


U T ' 
U T P 
U T P 2 


and DA P = A p P where D = in [65] , P is called the Stieltjes matrix 

associated to A P 


2.4.6 the production matrix 

In the context of ordinary Riordan arrays, the production matrix associated to a proper 
Riordan array takes on spécial form: 

Proposition 2.8 [2$j Let P be an infinité production matrix and let A p be the matrix 
induced by P. Then Ap is an (ordinary ) Riordan matrix if and only if P is of the form 


2.4. Riordan arrays and orthogonal polynomials 
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6 

Ct 0 

0 

0 

0 

0 •••" 


ai 

Cto 

0 

0 

0 ••• 

6 

et 2 

0\ 

Cto 

0 

0 ••• 

</» 

CO 

«3 

02 

0\ 

Cto 

0 ••• 

£4 

0^4 

a 3 

02 

0\ 

et 

£5 

«5 

Ct4 

Ct 3 

0-2 

Cti • • • 


Moreover, columns 0 and 1 of the matrix P are the Z— and A— sequences, respectively 
of the Riordan array Ap 


2.4. Riordan arrays and orthogonal polynomials 






Chapter 

Sequences of the orthogonal polynomials of 
Sobolev 

3.1 Définitions and properties 

We are given 2 f un otion al moments a and r and the bilinear form <f> defined by 

<i>(p,q) = (<*,p q) + (apV), p,q£V. ( 3 . 1 ) 

Here, we prove that, when 0 is regular, we can classify the polynomials solution of the 
Bochner’s differential équation. We point out that 0 is symmetric (4>(p,q) = 4>{q,p)) ■ 

Définition 3.1 [Tl] The numbers 

0n,m := 4>(x n , x m ), m, n > 0. 

are called the moments of 0 (.,.) , for any symmetric bilinear form 0(.,.). 

Définition 3.2 jTj. |37j A symmetric bilinear form 0 is said to be quasi-definite (resp 
defined positive) if it satisfies the Hamburger condition 

A n (0) =: det [0ij]" j=o ^ 0 (resp. A n (0) >0, n > 0). (3.2) 

Lemma 3.1 fÏÏÏ | / Let cj) be a symmetric bilinear form. For any normalized sequence of 
polynomials {\P n (x)} n>0 we hâve 

A„(0) = det {0(dh, ^)}A =0 ,n> 0. (3.3) 

Proof. For any n > 0, let A = (ajj)T= 0 be the square matrix of order (n + 1) and the 
sequence j’Fj(x) j defined by 

n 

= y ^a j>k ^k(x) , 0 < j < n. 

k =0 


3 
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Then, we hâve 


{#*,(*).»,(*))}• x* = U(Mx),*Ax)) 


and 


J t,j =0 


(x),^(x)) 


i,j=0 


i,j =0 


So, 


i , e 


<t>{*i(x), *j(x)) = A [<fi(Vi(x), %(x))]P =0 A 

J 1,3= 0 


det 


v ï ; j) . _ = (det A) 2 det [0(dk, ^j)] iJ=0 • 


jj=0 


If we choose a lower triangular matrix su ch that üjj = 1 and take ^(a;) = ad, for 
j = 0 , 1 , 2 , n, then 

A M = det [<j>{x\ ad )]" J=0 := det [ 0 ij]^ =o = det [ 0 (tfi, ^)]" j=0 . 

Hence, |3.3| holds. ■ 


Lemma 3.2 ^7 ]/ A symmetric bilinear form is quasi-defined (resp defined positive) if and 
only if there exist a sequence of polynomials {P n (a ;)} n>0 and real k n ^ 0 (resp k n > 0) 
for n > 0 such that 


0(P«, Pm) k n ô nm n , m ^ 0. 


(3.4) 


In this case, an SP {P„(a :)} n>0 
factor) determined. 


satisfying S.f is uniquely (to a non-zero multiplicative 


Corollary 3.1 The following statements are équivalent. 

a) {P n (x)} n> o is a sequence of orthogonal polynomials with respect to 0 

f 0(x m , P n (x )) = 0 ifO < m < n — 1, 

[ f)(x n , P n fx)) ± 0. 

b) For every polynomial i t(x) of degree m, 

j 0 ( 7 r(x), P n (x)) = 0 if 0 < m < n — 1, 

1 0 ( 7 r(x),P„(x)) 7 ^ 0 i/ m — n. 


We now investigate the regularity of the form 0 defined by 3.1 


Remark 3.1 The intégral représentation of a quasi-deûnite moment function is a.lways 
possible according to Duran’s classical theorem. It is convenient to dehne the orthogo- 
nality following in the sense of Sobolev. 


3.1. Définitions and properties 
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Définition 3.3 We say that {P n (x)} n>0 is a sequence of orthogonal polynomials of 

i.e., 


Sobolev (S O P S) if it is orthogonal to 0 defined by 3.1 


0(P n , P m ) := (cr, P n P m ) + (r, P' n P' m ) = kj n , m , with k n ^ 0, n, m > 0. 


(3.5) 


Lemma 3.3 JJ1\ \3l\ | Let 0 be a quasi-deftnite büinear form defined by 3.1 and let 
{P n ( x )}™=o be an SPOS with respect to 0 . Then, a is the canonical functional of 
{ p n{x)}n=oP-e-, 

(cr, 1) 7 ^ O.and (a, P n (x)) = 0, n > 1. 

We give necessary and sufficient conditions for the form 0 to be regular. To this end, 
we investigate three cases. 

Case 1: a is regnlar. 

Case 2: r is regnlar. 

Case 3: cr and r are not regular. 

Définition 3.4 The sequence of dérivatives of the sequence of polynomials {K n (x, y)} n>{] 
is defined by 


At ’* 1 (x, y) = 


d r+s 

dx r dy : 


K n (x, y), n> 0. 


Let us dénoté by{Q n (x)}, {R n (x)}^ L 0 and {Pn(x)}™ =0 corresponding S PO N veri- 
fying the orthogonality conditions 

(r, Qm(x)Qn(x)) = kn (t) ô m ,n 71,171 >0 

< (0, R m (x)R n (x)) = k n (0) n, m > 0 

(cr, P n (x)P m (x)) = k n (a) ô m , n n,m> 0. 

Now, we begin with the first case (cr is regular). where { S c ,p) := p(c), c G M 
Lemma 3.4 JJdj / The büinear and symmetric form 0 defined by 

0(p, q) ■= {°,pq) + a{àc,p'q'} P,qeV 

is quasi-definite if and only if 

1 + aK^-iic, c) 7 ^ 0, n > 1. (3.6) 


If 0 is quasi-definite then 


R n (x) = P n (x) 


aP' n fx)K ( ^(x,c) 
1 + aK^(c, c ) 


for n > 0 . 


(3.7) 


3.1. Définitions and properties 
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Moreover, with K_i(x,y) = 0 


or 


= 0 


k n (fi) = k n (cr) 

&TI— 1 

n > 0 

a n = 1 + aKn-i (c, c) 

n > 0 

*<*) = ™ 

1 

< 

n > 0 

P n (x) ^ R n {x ), 

\/n > 0 

k Poix) = Qo(x) et Pfix) 

= Qi(x). 


Now, we consider the second case (er = a<5 c ). 

Lemma 3.5 Ifityj Let fi be the bilinear and symmetric form defined by 


(j)(p,q) := a(S c ,pq) + (T,p'q '). 

Then, fi is quasi-defined (positive definite) if and only i/a / 0 and r is quasi-defined 
(resp. a > 0 and r is positive definite). 

In the case where fi is quasi-defined, we dénoté by{Q n (x)} and {i? n (x)}^ =0 respectively 
the corresponding SPON verifying the orthogonality 

{ (d~i Qm(x)Qn(x)') k n ( t ) 

{fii Rmix)Rnix)'} k-ri ( fi ) $m,n 

then 


Rn(c) 

= 0, n > 1 where ô c (x) is 

a canonical functional of {i? n (x)}^ =0 

fon-\-1 (^) 

= (n + l) 2 /c n (r), 

n > 0 


Qn{p^) 

= nh P ^ 

n > 0 


Rn(x ) 

= f c c nQ n _ 1 (x)dx, 

n > 1, 

Ro(x) = 1 


Proposition 3.1 Suppose that a^O and r is a moment classical functional verifying 
Pearson functional équation 

( h{x)T)' = h{x)r , 


3.1. Définitions and properties 
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or 


lj(x ) = y 'kjx 3 , i = 1,2. 
j =o 

Then {R n (x)}satisfies the non-homogeneous équation of order two 


with 


k(x)y" + (h(x) - k(x))y' + (nl'^x) - X n y(x)) = v n , 


A n = k 2 n(n - 1) + hrn 


Vn = k(c)Rn{c) + (hic) + l’ 2 (c)R! n (c ) for H > 0. 

3 Case cr, r are not quasi-definite and <f quasi-definite. 

Example 3.1 [3Ü] Consider the 2 functional moments a and r, defined by 

r = 


and (To = 3, a i = 1 et a n = -, n > 1, 

n + 1 

or {cr n } n>0 is the sequence of moments of FM a. 

As 


-, n = 0 


t(x — 


0, n > 1, 

we deduce that r is not quasi-definite. 

In addition, we hâve (cr) = 0, therefore a is not quasi-definite. 
Let us now deûne the bilinear form (j) by 


We know that 


< i>(p,q ) : = W,p q) + o(à<h'p'q r ), p, qeV. 


4>{p-,p) = / p 2 (x)dx + Gt/(0) 2 > 0; 


(3.8) 


then, 0 is positive dehnite =>■ <f> is quasi-dehned. 

We can construct more such examples via 

Lemma 3.6 [Jïtf For any symmetric bilinear form 'h, there is a positive-definite moment 
functional a such that the symmetric bilinear form 0 defined by 

0(p,g) := tf(p', q') + (a,pq) 


3.1. Définitions and properties 
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is positive-definite. 

In particular, for any FM t, there exists a positive- defined FM a such that 

4>(p,q) ■= q) + {r,pq) 


is positive-definite. 


3.2 Sobolev Orthogonal polynomials and Second or- 
der Differential Equations 


Remark 3.2 Let fi be a quasi-definite bilinear form given by the relation 3.1 where a 
and r are two FM. 

Our aim liere is to détermine necessary conditions for the SOPS with respect to fi, 
satisfy a second order differential équation having the form 

L [y] (x) = l 2 (x)y” + h(x)y' = X n y, with L [.] the differential operator (3.9) 
where the coefficients h, i — 1,2 and X n , n > 0 are as follow 


h(x) — E kjX j , * = 1,2 

j=o 


X n = n{n - l)l 2 2 + nln, n > 0 


k ^11 + ^22 0 . 


Theorem 3.1 [fiSj Consider the symmetric bilinear form fi defined by 3.1 The fol- 
lowing statements are équivalent. 

a) The differential operator L [.], defined by 3.9 satisfies the following relation 


fi(L[p\,q) = fi(p,L[q\), p,qeV. (3.10) 

b) The tuio functional moments a and r satisfy the following functional équations 

(l 2 (x)ay — l\(x)a = 0 (3.11) 

1 2 {x)t' — 1\{x)t — 0 (3.12) 

c) The moments of fi, a and r are given by 

fimn &m+n F U* ** T m +n—2, Tfl, H F 0 


(3.13) 


T-i = T_ 2 = 0, 


3.2. Sobolev Orthogonal polynomials and Second order Differential Equations 
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and 


(nl22 + hi)o'n+ 1 + { n hi + ho) <J n + nho&n-i — 0 , n > 0 , 
cr_ i = 0 


[{n + 2)Z 2 2 + ^ll] Tn +1 + [[n + 1)^21 + bo] T n + ^^20 r n-l — 0, U > 0 

T-i = 0 


(3.14) 


(3.15) 


Moreover, if<f>(.,.) is quasi-definite and {P n (x)}^l 0 is a SOPS, then the conditions 
a), b) and c) are équivalent to condition d), given below 


d) {P n (x)} n=0 satisfies the differential équation 3.9 


Définition 3.5 The functional équations 3.11|and|3.12 are the Sobolev weight équations 


of the operator L [.] with respect to the form </>(.,.). Equations 3.14 and 3.15 are the 
moments of Sobolev for L [.] with respect to 0(.,.). 


Remark 3.3 Theorem 3T is a generalization of the Krall theorem concerning the differ¬ 
ential équations of the 2 nd order of the sériés of orthogonal polynomials. 


Corollary 3.2 Let 0(.,.) be a symmetric bilinear form defined by 3.1 Then there exists a 


sequence of orthogonal polynomials of Sobolev with respect to </>(.,.) satisfying a differential 
équation of type \3.9 if and only if both of the following are satisfied 

1) </>(.,.) is quasi-definite 

2) The moments {o’ n }f =0 and {r n }^L 0 of a and r respectively satisfy the moment équa¬ 
tion 3.14 and 3.15 \ 


Remark 3.4 By deriving 3.9 with respect to x, and setting y'{x) = Z(x), one obtains 
the following second-order differential équation 


M [Z] (x) = l 2 {x)Z'\x) + (l' 2 (x) + ï{{x))Z\x) = (A n+1 - Vl(x))Z(x). 


(3.16) 


This leads us to dénoté |3T2] and |3.15| respectively by the weight and moment équations 
of M [Z] (x) . 


Proposition 3.2 Let L [.] be the differential expression defined by 3.9) and </>(.,.) is the 


bilinear form given in 3.1. If L{p) = A p and L{q ) = vq for each p,q eV and A, rel, 
with \ fi v, then 0(p, q) = 0. 

For any solutions a and r of the Sobolev weight équations 3.11 and 3.13\ respectively. 


Proof. This following from Theorem 3.1, one lias 

(A -v)(j>{p,q) = (j){\p, q) - 4>{p, vq) 


= <f>{L{p),q) - (j>(p,l(q )) = 0. 


3.2. Sobolev Orthogonal polynomials and Second order Differential Equations 
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i.e. 

00, q) = o. 


3.3 Classification of the Sequence Orthogonal Poly¬ 
nomial of Sobolev 


In this section, we consider the bilinear form fi defined by relation |3.1| Knowing that 0 
is a linear combination of two fnnctional moments, we will try to classify the sequences 


of polynomials that are solution of the differential équation |3.9[ For this purpose, it is 
assumed initially that one of the two FM is quasi-dehnite and in a second tirne that the 
two are not quasi-dehnite. 


3.3.1 a quasi-definite 


Theorem 3.2 


Suppose that there exists an orthogonal S PO S {P n (x)} n> 0 with 


respect to fi defined by 3.1 solution of équation 3.9. If a is quasi-definite, then {P n (x)} 


n> 0 


is necessarily classical with respect to a and t = kl 2 (x)a for some real constant k. 
In addition, we hâve r = 0 if k = 0, or else r is also quasi-definite. 


Remark 3.5 [49j If a is quasi-dehnite and r = cô a for ail c G M* and for ail a G M, 
then any OP S with respect to fi given by |3. 1| is called S PO S, but it does not satisfy the 


differential équation of Bochner’s form 3.9 


Moreover, this System satishes another differential équation of order 2 with the coeffi¬ 
cient dépendent on parameter n. 


Theorem 3.3 [J 7} / The classical orthogonal polynomials are also characterized by the fact 
that they are the only OPS satisfying the following condition: 

There exist 2 functional moments T\ not identically zéro and r 0 such that 

(ti, P'mK) + ( r 0, PmPn) = 0 , for m n and n = 0 , 1 . (3.17) 


3.3.2 r quasi-definite but a is not quasi-definite 

Theorem 3.4 J Suppose that there exists an orthogonal SPOS {P n (%)} n>0 with 
respect to fi defined by \3.1\ solution of equation \3.E\ If r is quasi-definite then 

i) m x)} n > i is a CO P S with respect to r vérifiés the differential équation 

ii) {P n (x )} n> o is a WOPS with respect to a. 

iii ) l 2 (x)a = kr, for some k E M. Furthermore, either l 2 (x)a = 0, or l 2 (x)a is 
quasi-definite. 


3.16 


3.3. Classification of the Sequence Orthogonal Polynomial of Sobolev 
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3.3.3 a and r are not quasi-definite 


Theorem 3.5 


Let {Pn{x )} n> o be a sequence of orthogonal polynomials with re¬ 


spect to (j) defined by \3.1\ If a and r are not quasi-defined, then the polynomial sequence 
{P n (x)} n>0 does not satisfy the Bochner équation. 


3.4 The d —Orthogonality of Sobolev 

We give a generalization of the orthogonality of Sobolev, starting with the study of the 
regularity of the bilinear forms and then giving the conditions on the polynomial sequences 
d—orthogonal in the sense of Sobolev. We also show that the classical sequences are not 
classical in the sense of the generalization of the Sobolev orthogonality. 

Proposition 3.3 Let {P n (x)} n>0 d—OPS with respect to regular U = (uo,u±, ...Ud- i) T ■ 
Let V = (vqjVi, ...u^-i) such as 


v r = u r + a r ô c , 0 < r < d — 1 

V regular if and only if 

1 + L*dn+k- l( c ; c ) 7^ 0 

where L n (x,c ) is the polynomial defined as 

k 


n—1 (d —1 

L dn+k(x, C) = Y 1 P ^ C ) Œ p p \ 

[ r=0 ^ Ur " 1 * d j+ r *3> 


OL r Pdj-\-r l , n / \ A Ql r Pin-\-r 


+ Pn(c) 


r _ q \0r , Pdn+r-Pn ) 

if a r = a then L* dn+k (x,c) = aL dn+k (x,c ) where L dn+k (x,c) defined in J73j 
Proof. Suppose that V is regular and let 


(3.18) 

(3.19) 

(3.20) 


m—1 


Qm Pm “1“ ^ ^ &m,iPi • &m,i ^ 


(3.21) 


i =0 


If we let m — dn + k with 0 < k < d — 1 , then we get 


("Uj-, Qdn+rPj) 


(U r , P d n + rPn ) = n” =o7 ° (î) _ 1)+r . +1 for j = 71 
OLdn+r,dj+r (u ri Pdj+rPj) for 0 if j 5; TL 1 


(3.22) 


and by |3.18 we also hâve 


{Uri Qdn+rPj ) 


(VriQdn+rPn) ^rQdn+ri.P)Pn^P) fo^ j H 
- aQdn+r{c)P n (c ) for 0 < j < n - 1 


(3.23) 


3.4. The d— Orthogonality of Sobolev 
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Accordingly, we obtain 


Qdn+k(%) Pdn+k{%) Q dn+k{.P) L d.n+k—1 (*^i c) 


Set 


x = c in 


3.24| to get 


(3.24) 


Qdn+k(c ) [1 + I/^ n+fe _ 1 (c; c)] = P dn+k (c) (3.25) 

where necessarily 1 + L^ n+fc _ 1 (c; c) ^ O,otherwise we also hâve P dn +k{c ) = 0. It now fol- 
lows by induction that c should be a common zéro for mor than d consecutive polynomials 
{P n } wich is impossible ( corollary 2.3 ) ■ 


Let us consider d— bilinear forms 0 O , 0i,..., fid-i (d > 1). 

Consider a symmetric d— dimensional vector bilinear forms 0 = (0 O , 0i, •••, 4>d-i) T 
defined by 

(j) r (p,q) := ( a r ,p q) + ( T ri pq ), ( 3 . 26 ) 

where a = (er 0 , cr 1; ..., a d - 1 ) and r = (r 0 , ..., T d - 1 ) are d—dimensional functionals. 

The d— orthogonality in the sense of Sobolev is a generalization of the orthogonality in 
the sense of Sobolev. On the basis of the d—orthogonality définition, we give the following 
définition. 


Définition 3.6 We say that the sequence {R n (x)} n>0 is a sequence of d—orthogonal 
polynomials of Sobolev (d — SOPS in short) if it is orthogonal to 0 which is defined by 


3.26 


i.e. 


4>a{Rrm Rn ) =0, Tl > md + « + 1, TU > 0, 
(0aj RmRmd+a ) 7^ d) ^7. A 0. 


(3.27) 


We give necessary and sufficient conditions for the vector form 0 defined by 3.26 to be 
regular. To this end, we study the following cases. 

case 1: a — (co, o'i, &d- 1 ) is regular. In this case, let {P n } n>0 be the d—COPS witli 

respect to a. 

We notice that if L n (x,y ) is an orthogonal polynomial for some n > 0, then its 
dérivative is given by 


n>0. 

Theorem 3.6 The symmetric d—dimensional vector bilinear form 0 = (0 O , 0i,..., (j> d -i) 
defined by 


(f) r {pi q) (cr r , pq) + a{5 c , p'q') p,q G V and 0 < r < d — 1 (3.28) 


3.4. The d— Orthogonality of Sobolev 
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is quasi- definite if and only if 


1 + aL^\(c, c) ^ 0, m> 1. 


(3.29) 


If 0 is quasi-definite then 


d aP m( x ) L dn+k-l(x,C ) fo Qn , 

RmyE) PmyE) /i i\ Tl ^ 0. (3.30) 

1 + aL) n _ ) l (c, c) 

Proof. For the first part , we assume that a is quasi definite. Then we can write 

m— 1 

Rm = Pm + a m,,iPi- (3.31) 


i =0 


If we assume that m = dn + k witli 0 < k < d — 1, then we get 
(<T r , R dn+r Pj) 


(ov, P dn+r P n ) for j = n 

® dn+r,dj+r ( T' ? Pdj+rPj ) for 0 T j T ^ 1 


(3.32) 


Frorn 3.28 , one lias 

(cr r , R dn + r Pj^ 

Accordingly, we obtain 

Rdn+k( x ) = Pdn+k{x) - aR' dn+k (c)c). 


(j) r (R dn+r , P n ) - aR' dn+r (c)Pf(c) for j = n 
-aR'dn + A c ) P n(c) for 0 < j < n - 1 


(3.33) 


(3.34) 


Differentiating |3.34 and setting x = c, then we get 

Rdn+k ( C ) 1 + p dn+k-l( C ’ C ) = P d n+ki C )i 

where necessarily l+L^+fc-, ( c ; c ) A 0, otherwise we also hâve P dn+k (c ) = 0. It now follows 
by induction that c should be a common zéro for more than d consecutive polynomials 
{P'n} m >î which is impossible. Hence, 3.30 follows immediately. ■ 

case 2: r = (r 0 , ri,..., r^-i) is regular. In this case, let {Q n } n>0 be the d— OPS with 
respect to to r. 

Now consider the following symmetric d— dimensional vector bilinear forrns 0 = 
(0o, 0i, •••, 0d-i) defined by 


0r(p, q) ■= ( a5 c ,p q) + ( T r ,pq'). 0 <r <d-l,ce 
Now we let {Qn(x)} n>0 be d-OPS with respect to r satisfying 


(3.35) 


vr, QnQm ) 


0 if n > dm + r + 1 

k n {r r ) if n = dm + r, m > 0 
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and {R n (x )} n> o be d- OPS with respect to 0 satisfying 


Vn RnRrri ) 


0 if n > dm + r + 1 
k n ((f> r ) if n = dm + r ,m > 0 


The next proposition will prove useful later on. 


Proposition 3.4 The vectorf is quasi-definite if and only if a ^ 0 and r is quasi definite. 
Then, 

R n (c) = 0 for n > r + 1, 
k n +i{4>) = {dm + r )mk dm+ \(r ), 

Qn = 

Proof. Assume that <f> is quasi definite 


0r {Rmi Rr 


On the other hand 


0 if n > dm + r + 1 ,m > 0 
k n {4>) if n = dm + r ,m > 0. 


(j) r ( R ,n ■ R n 


{aS{x c),R m Rn) T {r r iR ln Rn') 

^•Rm{^)Rn{c) T (j~ri R'm^'n') 


In particular, for m — 0 we hâve 


0r ( f M ). A/, ) aR n (c) 


0 if n > r + 1 
k n {4> r ) if n = r 


z.e., 


kr{4>r) ^ n 

a = — ; , 7^ 0 


R r {c] 

Rnfc) = 0 
Thus for m and n > 1 we hâve 


for 


n > r + 1 


(p r (R r n i Rn) {Tri R-m ) 


0 if n > dm + r + 1 ,m > 0 
k n ((f) r ) if n = dm + r ;m > 0 


Then, r is quasi definite and {R' n (x)} is d—OPS relative to r. In this case we get 


R' n {x) = nQ n - 1 , n > 1 


and 


(x r , R m R n ) (r r mnQ m ~\Q n—i) 


0 if n > dm + r + 1 

{dm + r)mkdm+r{ T r) if n = dm + r, m > 0 
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Now, we assume that a ^ 0, r is quasi definite and {H n } n>0 is MPS, such that H' n = 
nQn-i- This leads to 


MH m , H n ) = aH m (c) H n (c) +{r r , H' m H’ n ) 


a if m — n — 0 
0 if 77i = 0 and n > r + 1 

< 

0 if n > dm + r + 1 

(dm + r)mk dm+r _i(r r ) if n = dm + r ,m > 0 


i.e. {H n } n>0 is a d- OPS with respect to (f> and by nnicity, we obtain H n = R n . This 
complétés the proof. 

In this cas, we shall prove that the 2-classical orthogonal polynomials with respect to 
r, satisfy homogeneous third order lincar differential équation of Bonkhemis type. 


a n,i( x )Q ( 'n+z( X ) + a n,3( X )Qnl3( X ) + a n,2(x)Qnl 3 (x) + 07 n ^(x)Q n+3 (x) = 0, (3.36) 

with a Ui i(x) for i — 1, ...,4 are polynomials of degree less than or eqnal to i and ■ 


a n , 4 X ) = F nA (x)S 3 (x), 

is not define a 2—classical OPS in the sense of Sobolev 


Theorem 3.7 Assume that a^O and r is 2-classic moment functional vector. 

If dega ni j(x) < i — 1 then the sequence {.R n (a;)} n>0 satisfîes the following non- 
homogeneous third order linear differential équation 


a n- 1,4^1+3 + ( a n-l,3 “ a n-l,d^n+3 + K- 1,4 + a n-l,2 ~ Q^- 1 , 3 )-^1+3 + 
( a n- 1,1 - «n— 1,2 “ a n- 1,4 + a n-l, 3 )-^n+3 = v n, 


(3.37) 


where 


V n — Q!ra—1,4-R n -|_3+ (û^n-1,3 a n-l,4)-^n+3“i" ( a n -l ,4 + a n-l,2 a n-l,3)-^n+3\a:=c, n — 0 (3.38) 

If dega n j(x) = i then the sequence {iî ri (m)} ri>0 satisfies the following homogeneous 
four order linear differential équation: 

«n—l,4-Ri+3 + «n-1,3^1+3 + «ra-1,2^1+3 + a n_l,l ^,.+3 = 0- (3.39) 


Proof. Let {Q n } n>0 be a classical orthogonal polynomial sequence with respect to r 
satisfying 3.36 and Rl n (x ) = nQ n _ i, n > 1. 


For n = 0, i?o satisfies 3.37 
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For n > 1 and dega ni i(a;) < i — 1 one has 


(«ra—l,4-Ri+3 + (oin- 1,3 “ a ' n _l, 4)^+3 + («n-1,4 + a n-l,2 ~ «n-l^-^+s)' 

= (n + 4)(o; n — l,4<3i+2 + («n-1,3 ~ a îi-l,4)Qn+2 + («ra-1,2 + «„-l,4 — 

= (n + 4) (a„_i, 4:Qn+2 + a n-l,3<3n+2 + «n-1,2<5n| 2 + («n-1,2 + «1-1,4 _ «n-l, 3 )^+ 2 ) 

= (n + 4)(—Oin-i^Qn-l + («n_1,2 + «1-1,4 — «n-1, 3 )^+ 2 ) 

= + a^_i i2 + ai 3 2i i4 - Vl,3)^n+3 


?( 2 ) 


>( 1 ) 


and 


( 

«n-l,4-R n+3 + K-1,3 - a îi-l,4)-^n+3 + («n -1,4 + a n- 1,2 ^n-l,3)-^n+3 

= f(—a n - 1.1 + «1-1,2 + «1-1,4 _ a n-l,3)K+3^ 

= (— «n-1,1 + a ^_ i j2 + « 1 - 1,4 _ a n-l, 3)^+3 + V n 


?( 2 ) 


( 1 ) 


î.e. 


«n-l,4-Rl+3 + («n-1,3 ~ «n-1,4)-^1+3 + («n-1,4 + «Ji-1,2 ~ «n-1, 3 )-^l+3 + 

/ „,(3) 


>( 2 ) 


,(1) 


(«n—1,1 


a 


n—1,2 


a, 


1 -1,4 4“ «n-l,3)-^™+3 v r 


for some constant v n . Since R n (c) = 0, n > 1 we hâve 3.38 


If dega n ,j(a:) = i and according to 3.40 we obtain 3.39 
Tlins the proof is complété. ■ 


(3.40) 
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Characterizations of classical d-OPS 


4.1 Introduction 

Polynomials are spécial functions that we can differentiate and integrate many tirnes. 
Any polynomial may be written as a linear combination of ail other polynomials of lower 
degree, i.e., each polynomial of degree n mnst satisfy a linear récurrence of order n (n + 1 
ternis). If the récurrence is Unité, i.e., there exists m > 0 such that 

n— 2 

Pn (*^) (d-7 iX b n ) P n —\ H - ^ ^ Q j n,kPki 

k=m 

we say that the sequence of polynomials {P n } n veriûes an orthogonality of such type well 
dehned. 

Ail notions of generalization of the orthogonality introduced by the authors, i.e., mul- 
tivariate orthogonality, multiple orthogonality, bi-orthogonality, d-orthogonality, ..., are 
obviously for searching a solution of some problems like Padé approximation, rationality 
or irrationality solution of higher order ordinary or partial dif férent!al équations, ... 

Classical orthogonal polynomials play an important rôle among orthogonal polynomi- 
als. The concept of d-orthogonality, d being a positive integer, appears as a particular 
case of the general multiple orthogonality. In particular, they satisfy orthogonality con¬ 
ditions with respect to several positive measures. Some recent works were focused on the 
analysis of properties of d-orthogonal polynomials generalizing the classical orthogonal 
polynomials and some characterization theorems were derived. 

Many examples of multiple orthogonal polynomials as well as of d-orthogonal polyno¬ 
mials were derived. They satisfy a linear differential équation of order d + 1, i.e., when 
d = 1 we meet the standard orthogonality. In this case, we know that there exist only 
four families of continuons and four families of discrète classical orthogonal polynomials 
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solution of the Bochner’s équation, i.e., a linear équation of second order of the form 


a 2 (x) y" + pi (x) y' + y = 0, 


(4.1) 


where «2 and Pi are polynomials of degree two and one respectively and where 7 is 
constant. 


The solution of the équation (4.1 ) are known as classical orthogonal polynomials. Start- 
ing frorn tliis point of view, many characterizations of the classical orthogonal polynomials 
are given. 


In this chapter, we wish to generalize most of these well-known characterizations for 
standard orthogonality in the case of d— orthogonality. Most of the results are given in 
the case d, but according to the complexity of the calculations we give some results in the 
case d = 2. 


Notice that the application of the classical polynomials is very large and in various 
domain. This is our wishes for this study. We hope that our results enlighten such a 
way for the applications in a new domain that not possible in the case d = 1, at least at 
the moment, we know that the classical 2 -orthogonal polynomials satisfy a third order 
differential équation whicli is impossible by the standard orthogonality. 


4.2 Preliminaries and notations 

4.2.1 m-Symmetric sequences 

Problems related to symmetrization of sequences of orthogonal polynomials on the real 
line play an important rôle |24j . In general, we hâve 

Définition 4.1 [füj Let mbea nonnegative integer. 

(a) A polynomials sequence {P n } n>0 is called m-symmetric if 

P n (u>k'x) = w™P n (x) , n > 0, where uj k = exp {2k-Ki/ (m + 1)} , 1 < k < m. 

(b) A form T = (To,..., Td-i) is called m-symmetric if 

(T t; ,x (m+1)n+#1 ) = 0, 0 < y < m, v ^ y, n > 0, 
for each 0 < v < d — 1. In this case, 

(T^, P (x)) = ui% (un, P (idkx)) , where uj k — exp {2km/ (m + 1)} . 
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For the particular case: m — 1, we meet the well known notion of symmetric PS EU- 
The components of the sequence {P n } n>0 are (m + 1 ) sequences {P%} n>0 , 0 < n < m, 
such that 

P(rn+\)n+,i. = X^P% ( X m+1 ) , 0 < fi < m, U >0. (4.2) 

Définition 4.2 [29j A d-OPS {P n } n > 0 is called classical according to the Hahn property 
if Qn — (n + l )” 1 P' n+ 1 is also a d-OPS. 

Since the dérivative operator takes each d-symmetric PS into a d-symmetric PS, then 
a d-OPS {P n } n >o is called d-symmetric classical, if Q n is d-symmetric d-OPS. 

A characteristic property of m-symmetric PS is given by the following 

Proposition 4.1 f29. j, [16] For each d-OPS {P n } n>0 with respect to U, the following state- 
ments are équivalent. 

(i) The formU is d-symmetric, 

(ii) The sequence {P n } n>0 is d-symmetric, 

(iii) The sequence {P n } n>0 satisfying the récurrence 

Pn+d+ î {x) = xP n+d (x) - 7 ° + i-P n {x ), n > 0 , 

(4.3) 

P n (x) = x n , 0 < n < d. 

Note here that if {P n } n>0 a m—symmetric d—OPS, ail the components {Pf} n>0 , 
0 < n < m, are d—orthogonal and if moreover {P n } n > o i s classical, then ail the (d + 1 ) 
components are classical. 

4.2.2 Characterization of m— symmetric 

Let d be a positive integer and m de a nonnegative integer satisfying m < d. Next, 
we give a necessary condition on m and d to hâve an m— symmetric d—OPS and two 
characterizations of m— symmetric d—OPS. We dénoté by X k the multiplication operator 
by X k in V 

Theorem 4.1 flûjLet {P n } n>0 be a d— OPS . Then the following properties are équivalent: 

1. The PS {Pn } n > 0 is m—symmetric 

2 . d + 1 is a multiple of m + 1 , say d + 1 = p(m + 1 ), and the PS {P n } n>0 satisfies a 
(d + 1 )—order récurrence relation of type 

p 

XP n = P n +1 + 7 n,jPn- j(m+ 1)+1 

3 = 1 

with 7 njP ^4 0 and the convention P_ n = 0 for ail n > 1 
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Theorem 4.2 [îW Let {P n } n> o be a d—OPS . Then its components {P[} n>() k — 0, 
are d— orthogonal 

Classical d— OPS 

Theorem 4.3 [TE] If {P n } n>0 «s a d— symmetric classical d—OPS, then its the compo¬ 
nents {Pn} n>0 k — 0, ...,d, are classical d— orthogonal 

Theorem 4.4 [TB] Let {-Pn} n>0 be an m—symmetric classical d—OPS then its first com- 
ponent {Pn} n>0 is a classical d— OPS 


4.2.3 d-quasi-orthogonality 

In the rest of this section, we recall the notion of d-quasi-orthogonality and some it 
characterizations. 


Définition 4.3 [54] A sequence {P n } n>0 is said d-quasi-orthogonal of order s with respect 
to T = (r 0 ,r rf _!) T , if for every 0 < a < d — 1, there exist s a > 0 and cr a > s a integers 
such that 


(T a ,P m P n ) = 0, n > (m + s a ) d + a + 1, m> 0, 
(To) P<r a P(cr a +s a )d+a ) 7 ^ 0, 7TI ^ 0, 


(4.4) 


for every 0 < a < d — 1. We put s = maxs Q , 0 < a < d — 1. 


Définition 4.4 [54]A sequence {P n } n>0 is said strictly d— quasi orthogonal of order s 

T 

with respect to T = (r 0 ,..., r^_i) if it satisfies 


(T v ,P m P n ) = 0, n > (m + s)d + v + 1, m > 0, 

(T a ,P m P {m+s)d+v ) 0, m 0, 
for every 0 < v < d — 1 with s = maxo^^d-i s a 


(4.5) 


Définition 4.5 [5l] Let v G V and d,s > 1 integers. The sequence {-P n } n>0 is strictly 
s/d-orthogonal with respect to v if it fulfills: 


(v , PmPn) =0, n > md + s, m > 0, 

(fo PmPmd+s—l ) 7^ Oj TTl 0, 

With this définition we can express: {Pn} n> 0 is d-OPS with respect to 

(fo,..., Pd-i) T if and only if is strictly (a + 1) /d-orthogonal with respect to r Q , 


(4.6) 

r = 
o < 


a < d — 1. 

A d -OPS {P n } n >o with respect to fo,can be d—quasi orthogonal of order s with respect 
to — of dimension d . We hâve the following resuit. 
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Proposition 4.2 [Tfij Let {Pn} n>0 a d -OPS with respect V . Suppose that exist T of 
dimension d such that {P n } n>0 is quasi orthogonal of ordre s with respect to T. Then 
necessarily d > d 

We consider only the case d = d , where we need the following characterization . 

Theorem 4.5 [5fÿ For any sequence {P n } n>0 a d -OPS with respect V , the following 
statements are équivalent 

(a) There exist T of dimension d such that {P n } n>0 is quasi -orthogonal of order s with 

respect to T 

(b) There exist T of dimension d and s a , cr a > s a integers, 0 < a < d — 1 such that 



(4.7) 


for each 0 < a < d — 1 


(c) There exist — of dimension d and d 2 polynomials (p'f,0<u,a<d—l 


d -1 



(4.8) 


where 


deg 0" = s a , 0 < a < d — 1, and if d > 2, 

deg^“ < s a , 0 < a < u — 1, if l<a<d—1 

deg </>“ < g — 1, u+l<u<d — 1 if 0 < a < d — 2, 


(4.9) 


(d) There exist T of dimension d such that is strictly quasi orthogonal of order s with 
respect to T = (T 0 , ...,r d _i) T 


(e) There exist —of dimension d and s a > 0 ^integers, 0 < a < d — 1 such that 


(r a , P n ) = 0, n > s a d + o + l, 
(r ai P Sa d+ a ) -j- 0, 


(4.10) 


for each 0 < a < d — 1 
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Lemma 4.1 J A sequence {P n } n>0 satisfy 4-4 can n °t be d— quasi-orthogonal of order 
s — 1 with respect to T = (T 0 ,r rf _!) T 


Définition 4.6 [20j A sequence {P n } n>0 is said to be weakly d-orthogonal of index ( k , l) 
with respect to Y = (r 0 , ...,rd_i ) 7 if there exist integers, k,l > 1 such that 


(r a , P n ) = 0, n > d (k a - 1) + a + 1, 
(r a , Pd{k a -l)+a) ~f~ o, 


(4.11) 


j (T a ,xP n ) = 0, n> d(l a -l) + a + l, 

y {r a ,xP d{ i a _ 1)+a ) 7 ^ o, 

with k = max k a and l = ma xl a for every 0 < a < d — 1 . 


(4.12) 


Then, a strictly d-quasi-orthogonal sequence of order s with respect to T, is weakly 
d-orthogonal of index (s + 1, s + 2) with respect to U. 


Theorem 4.6 For each sequence (P n } n>0 d-OPS with respect to —, then the follow- 
ing statements are équivalent. 

(i) There exists C G V and an integer s > 1 such that 

(C,P n ) = 0, n>s and (£, P s _i)^0, (4-13) 


(il) There exist C G V, an integer s > 1 and d polynomials <f a , 0 < a < d — 1 such that 
C = <f> a u a with the following properties: if s — 1 = qd + r, 0 < r < d — 1, we 

hâve 

deg(p r = q, 0 < r < d — 1, and if d> 2, 

deg0“ < q, 0 < a < r — 1, «/ 1 < r < d — 1, (4-14) 

deg0“ < q — 1, r + l<o<d — 1, if 0 < r < d — 2, 

Another characterization of d-quasi orthogonality in terms of sequences is the following 


Proposition 4.3 J TÏÏjFor any two d—OPS,s \P n } n> 0 an d {Qn} n>0 relative to T and V 
respectively , the following are équivalent 


a T and V satisfy the relations 2.15 and 2.16\ for some non zéro d x d matrix polyno¬ 
mials © = 


b there is nonnegative integer l such that 


dl 

Pn{x) = Qn(x ) + y; a n)i Q n -i where a n>d i ^ 0 (4.15) 

i =1 
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4.3 Classical <i-orthogonal polynomials 

The characterizations of classical orthogonal polynomials in the algebraic aspect made by 
Maroni in 1533. can be broaden here in the sens of d-orthogonality. We point out that 
the resnlts given in the theorem below were the first characterization of d-classical OPS 
pioneered by Maroni and Douak (see, e.g., [30]). Here we first aim to simplify their proof 
and we shall give another characterization of d— orthogonality. 

Theorem 4.7 For any d-OPS {P„} n>0 with respect to T = (r 0 , ...,Pd_i) 7 , the next 
statements are équivalent: 


(i) The d-OPS {P n } n>0 possesses the Hahn’s property. 

(ii) There exist two d x d polynomials matrix $ and T such that 

($r)' + Tr = 0, 

where $ = (0))) and T are in the following forms [30] 


T = 


0 

0 0) fi 


0 

0 

0 d- 1 

Ù -1 


with 0 (x) is a polynomial of degree 1, {0u})!=i are constants and 


$ = 


0o fa) 

0? 0) 


r 1 (*) 

^ (*) 


'd -1 


(X 


id-l 

J d-1 


X 


where (f) v a are polynomials such that: 

deg 0" < 1, 0 < v < a T 1, if 0 < a < d — 2, 

deg f>a = 0, a + 2<v<d — 1, if 0 < a < d — 3, 
deg , < 2 and deg0)]_ 1 <1, 1 < v < d — 1. 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


(iii) {P n } n> o is d-quasi-orthogonal of order 2 at most with respect to V = 4>r and {Q n } 
is d-OPS with respect to V = $r. 
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Proof. we prove the implication (i) =>• (ii) . Suppose that {P n } n>0 is classical d-OPS 
and {Q n } d-OPS with respect to V = (v 0 , r>i, ..., i) T where Q n = (n + l)^ 1 P' 1+1 - By 
theorem (|2 .11) we hâve 


V dn+n = fin^Vo + i + ... + i 


then, taking the dérivative and using again the System (2.16) and using also (2.3) we 
obtain 

- (dn + /i + 1) T dn+fl+1 = v 0 + ui + ... + {(f) d n -yy v d -1 

-dK^To - X ):iZo [(* + !)</C + <%i+ i<?] r m, 

Fdn+fi+l = Vn,!i+ 1^0 + Vn,n+ 1^1 + ••• + ^nj+l^d-l ■= ^+ 1^0 + + ■■■ + ^+ 1^-1 

(4.20) 

which implies 

Vq + (<Pn tli yv 1 + ... + {tâ^y V d _ 1 = [d^V - (dn + /i + 1) pj +1 ] u 0 + 

Ef=o [(* + !)<„ + dÇi+itâJ - (dn + fi + 1) 

note by 

k 0 (n, /i) : = d(j ) d “ V ~ ( dn + V + 1) ¥>J+i 
and 

fca(w, h) : = OC" 1 + d^ a (f) d ~y -(dn + fi + 1) <p* +1 

with 


degfc 0 (n,/i) = 


and 


< n if 0 < /i < d — 2 

< n + 1 if fi — d — 1 


<n — 1 if /r + 2<a<d — 1 and 0 < fi < d — 3 
deg [fc a (ra, /i)] = { < n if 1 < a < yU + 1 and 0 < p, < d — 2 
= n if /i = d — 1 and 1 < a < d — 1 

If we take n = 1, n — 0 and (4>{j) = 1, we hâve 

d -1 

«o = fcn(1.0)«. + y^Ul,QK. 


(4.21) 


a=l 
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And for n = 1, /x = 1, we obtain 

d— 1 

(0?,i)'mo + ^1 = ^o(l, l)wo + X] M 1 ; i)^-, (4.22) 

c*:=l 

to simplify calculations take (01 J ) / = 1 for i < j 

i.e. 

vi = [k Q ( 1,1) - fc 0 (l, 0)] u 0 + E [fc Q (l, 1) - k a ( 1, 0)] u a 

q=i (4.23) 


put n = 1, 1 < /i < d — 1, we obtain 


d -1 

v^= [k 0 (l,n)-k 0 (l,n-l))u 0 +^2[k a (l,n)-k a (l,n-l)]u a ( 4 . 24 ) 

Oi=l 

Note by 

k a { 1,0):= 0o a := 0,..., d — 1 

[k a ( 1, /i) — k Q ( 1, /x — 1)] := 0" a := 0,..., d — 1 and fi = 1, ...d — 1 

with 


deg 0" < 1 , if 0 < a < ytx + 1 and 0 < /x < d — 2 
deg 0“ = 0 if / u + 2<o<d—1 and 0 < /x < d — 3 
deg0^_ 1 < 1 if 1 < a < d — 1 and deg0°_ x < 2 
hence we hâve obtain 

v = 4>(x)r. 


No, by taking n — 0, /x 


d — 1 in (4.20) and using also (2.3) we obtain 


— dT0 — CoIA + CiT2 + ... + Cd-2^ d-l + Cd-lV , d _ 1 

— —d [0iTo + AiTi + À 2 T 2 + ... + Àd-iT^-i], deg0i = 1 


i.e 


v 


t 

d— 1 


-[—d0iTo — (co + dAi)Ti — .... — ( Cd -2 + dA<2_i)rk_i] 

Cd-l 

= — [0r o + £ir 1 +.... —hCd-iTti-i] 


whence (4.16) and (4.17). 
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(ii) =^(iii) Since V = 4T, then 

(Va, Pn) = Z) ( r M> =0, Tl > d + fl + 1, 

fi=0 

(v a ,p d+a ) = (r a ,(j)^p d+a ) ± o, 

for 0 < a < d — 2 and 

(n d _i, P n ) = <r 0 , + EÎ=Î ( r /^ <f> d -iPn) = 0, n > 2d + n + 1, 

(v d -l,P 2 d) = (r 0 , </>°_i-P 2 d+/i) 7^ 0 . 

Then there exists 0 < t < 2 such that {P n } is d-quasi-orthogonal of order 2 — t with 
respect to V = $—. 

(iii) ^(i) Trivial since by hypothesis both {P n } n>0 and {Q n } n>0 are d-OPS. ■ 

Unlike in the case d = 1, if {P n } is d-classical (d > 2), we can’t conclude from Pearson 

équation that {Q n } is also d-classical. Indeed, since {Q n } is d-OPS with respect to 
V = $r, then 


($V)' = $'V + «hV 7 = $'V + $Tr := TiV 4» 

Then, the following problem is open: the dérivative of classical d-OPS is again classical 
d-OPS? 

We try to give an answer to this question in the présent section. 


Fairly, it seems as the first look that the answer is négative! In fact, we hâve an 
example p.86, Tableau II]. In spite of that the family (E) satisfies a second order 
linear differential équation, meanwhile this family isn’t classical at ail. 

Our results show that the solution of Boukhemis and Zerouki differential équation P3] 
are the only 2 -classical OPS for which any dérivative of that sequence is again 2—OP S. 
For this end, we set the following définition 


Définition 4.7 A sequence of d-orthogonal polynomials {P n } is called very classical d- 
OPS if {P n } as wcll as it’s dérivative of any order are classical d-OPS. 


So now, the question is: déterminé ail the very classical d—OPS. 

On the other hand, note here also that we can extract a simple characterization from 
proposition 2.6 as follows (see m) 


Proposition 4.4 Let {P n } n >o 


monic sequence verifying (4-4) su °h that 


(a) Xn+d-i,v = 0, 0 < v < n - 1 and Xv+d-i,v 7^ 0, v > 0, 


(h) {Qn } n> o ts d-OPS with respect to V. 
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Then, the sequence {P n } n>0 is d-classical OPS. 


According to the last assertion of Proposition 4.3, since {P n } n>0 is d-classical and 


Qn = {n + 1) _1 P n+ i is also d-OPS, and from theorem 
quasi-orthogonal of order two at most, then we hâve the following characterization of 
d-classical OPS 


4.7 


n> 0 

we know that {P n } n>0 is d- 


Corollary 4.1 Let {P n } n>0 be d-OPS with respect to T = (r 0 , ...,Pd_i) T . Then {P n } n>0 
is d-classical if and only if there exists an integer p < 2d such that 

n 

P n (x) = ^2 ^n,vQu (x ), n> p. (4.25) 

v=n—p 

Note that when d — 1, this characterization is given in [521 [53]. And for d — 2, this 
resuit proved ûrst by Maroni in [[58] using another approach. 

The next characterization is presented for the semi-classical case in mi, we give here 
a direct proof for the classical character and it can be also used to prove a similar resuit 
in the semi-classical case|71j . 


Theorem 4.8 f7I{ / Let {P n } n>0 be d-OPS with respect to T = (r 0 ,..., Pd_i) T . Then 
{Pn } n> o d-classical if and only if {Q n } n>0 is weakly d-orthogonal of index (1,2). 

Proof. Suppose {P n } n >o d-classical, then {Q n } n>0 is d-OPS with respect to V = <h—. 
Afterwards, we hâve 

(v^, x m Q n ) = 0, n > dm + p + 1, 


(v li ,x m Q dm+li ) J- 0. 

In particular, respectively for m = 0 and m — 1, we hâve 

{v^, Q n ) = 0, n > p + 1, 

{v», QQ ~f~ 0, 

and from the ( d + 2) order récurrence of {Q n }, we hâve 

(vp, xQ n ) =0, n>d + p + 1, 

(v^xQ^) = 7° +1 (v^QQ ± 0 . 


Comparing with the Définition 4.6, we conclude that (k,l) = (1,2). 

Conversely, suppose that {Q n } is weakly d-orthogonal of index (1,2) with respect 
to V and we show that {P n } is d-classical with respect to T, i.e., we will prove that 
T satisfy Pearson équation (4.16). Let W = (w 0 , ..., Wd-\) T such that W = V', i.e. 
(w a ,ir) = (v' a ,n) , V7T G V. Then 

{w a ,P n ) = -n{v a ,Q n -i) = 0, n>a + 2, 


(w a , P a+ 1 ) = -(« + !) (v a , Q a ) ^ 0. 
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Then, there exists a < t a < a + 1 such that 


Hence, by theorem 4.6 


(w a ,P n )= 0, n>t a + 1, 

(w a ,P ta ) ~f~ 0. 

there exist d polynomials a < fi < d — 1, such that 


w 


a 



with, if t Q = q a d + r Q , 0 < r Q < d — 1, we hâve 


deg'i/j^* = q a , 0 < r a < d — 1, and if d > 2, 

deg ^<q a , 0 < n < r a - 1, if 1 < r a < d - 1, 

deg î[}£ < q a — 1, r a +l<n<d — 1, if 0 < r a < d — 2. 

Since a < t a < a + 1, we distinguish two cases: 


(1) If 0 < a < d — 2, whence a<t Q ,<a + l<d. Consequently, q a < 0. 

(2) If a = d — 1, then d — 1 < td-i < d, whence qd-i < 1. 


If = 1, then r d _! = 0 and td-i = d, and if qd-i = 0, then r d -\ < d — 1. Wliich 
prove that ail polynomials are constant except t/’d-i = ax + b. 

Now, it remain to show that the matrix 'h is of the form (4.17). Indeed, from (2.3), we 
hâve 


v = 


1 

^ o 


v[ 

= — 

1 

7 

^ "C3 


), we obtain 


1 

o 

0 


0 

2 


0 

0 

d 


Ti 

r 2 


r d = </>?,() O) «0 + 5^ 01,0 ( x ) r /i 

z — * fl=l 

where 05,o ( x ) = ex + d and where (x) — £ for 1 < n < d — 1. Hence 


a—i 


0 1 0 ■■■ 0 

002 : 

; o 

o • ■ • 0 d-1 

0 (x) Ci 6 • • • Cd -1 


r 0 

Ti 

r d _i 
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Now, we let the form W = (wq, ...,Wd-i) T such that W = (xV)' , i.e. (w a ,7r) = 
{(xv a )' , 7r) , W G V. So, since {Q n } is weakly d-orthogonal of index (1,2) with respect 
to V, then 

(w a , P n ) = (( xv a )' , P n ) = —n (v a , xQ n _ i) = 0, n > d + a + 2, 
k, ^Pd+cH-i) = - (d + a + 1) (u a , xQd+c) 7^ 0. 

Then, there exists a < t a < d + a + 1 such that 

(w a ,P n )=0, n > + 1, 

{w a ,P ta ) 0, 

In addition, there exist d polynomials a < /x < d — 1, such that 



put 1 Q = q a d + r a , 0 < r a < d — 1, we hâve 

deg = g a , 0 < r a < d — 1, and if d > 2, 
deg < q a , 0 < n < r a — 1, if 1 < Tq, < d — 1 

deg < ç a — 1 , r a + l</i<d — 1 , if 0 < r a < d i— 2. 

Since a<t a <d + a + l, we distinguish several cases: 

(1) If 0 < a < d — 2, whence a < t a < d + a + 1 < 2d. Consequently, q a < 1. 

(2) If q a = 1, then t a = d + r a < d + a + 1 hence Tq, < a + 1. If = 0, then necessary 
a < r a . 

(3) If a — d — 1, then d — 1 < td-i < 2d, whence < 2. 

If = 2, then r^-i = 0, and if rd -1 < d — 1, then q a < 1. 

If = 0, then r^-i = d — 1. This show that only d°_ 1 is of degree two. 

Now, we hâve V = TT and ( xV)' = 0T. If we let V = $r, then 

0 - = (xV)' = xV + V = xd/T + $T. 

Hence $ = 0 — xT with V = and (d’h) 7 = TT. ■ 

Now for any sequence of OP S it’s necessary to know thcir explicit expression of the 
measure for which the orthogonality is satisfied. It is however good enough to give a such 
tool for determining the orthogonalization measure. The quasi-monomiality principle is 
a very useful tool as well. Because it’s not required any information about the classical 
character. Our approach deals only with the classical OPS and it will be possible to 
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generalized for semi classical case. In the following, we prove a new characterization of 
2—classical OPS, which will be without a doubt the pillar in this context. Note here that 
the next results are used before by Maroni and Douak to give an intégral représentation 
for sonie particular class studied starting from the corresponcling generating functions. 

Theorem 4.9 For any 2 -OPS {P n } n > 0 w ^th respect to T = (To, Ti) , the next statements 
are équivalent: 


(i) The functional momentT satisfy the Pearson équation (4-16) with (f.ll), (i.e. T is 
classical), 


3 3 

(ii) There exist {a l n } i=zl and {(T n }, J=x polynomials of degree n, and /i i; H 2 , K 2 , M 2 and 
N 3 ,L 3 of degree <1,2,3 respectively, such that 


with 


and 


al (x) r" + a\ (x) r' 0 + al (x) r 0 = 0, 

Pl (*) r'/ + PI {x) r; + p\ {x) r, = o, 


(4.26) 


h 2 (x) t 0 = k 2 (x)t 1 + l 3 (x) r;, 
(x) r 1 = m 2 {x) r 0 + n 3 ( x ) Tq. 


(4.27) 


Ho = u'<f) - u (0' - if), K 2 -- 
/il = 0lTi — Tt9\, M 2 = TT (<// 


Proof. Write the matrices in (4.16) when d 


1X9 1 - (j) 7Ti, L 3 
-if)-9J, N 3 

2 in the form 


= 0j9 — (j)TT, 
= 7 T(j) — w9, 



0 1 

and $ = 

U! ( X ) TT ( X ) 

y = 


if(x) £ 




with £ is constant and the polynomials i/j, ui, tt and 9 are of degree 1 at most and deg <f < 2. 
In this case, the System (4.16) is équivalent to 


o/To + aPo tt i -\- ttY^ — 0, 

((p - v») r 0 + PK + 0 ii\ + 9VI = o, 


(4.28) 


with tti = H - 1 = C te and 0 1 = 0 , -£ = C 


ite 


First we establish ail possible cases of the second équation in (4.26). So, eliminate Tg 

h 2 (x) r 0 = k 2 fx) Tl + l 3 (x) r;, ( 4 . 29 ) 


from (4.28) we get 
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where 


H 2 = ü/0 — UJ (0' — 0) , K 2 = u 6 i — 07Ti, L 3 = U)9 — 07T, 


(4.30) 


with deg II> < 2 , deg A4 < 2 and deg L 3 < 3. 

We distinguish three possible cases according to the degree of the polynomial H 2 (x). 
1) degH 2 = 0, i.e. H 2 (x) = C te := C. In this case, replacing r 0 and Tg obtained from 
(4.29) in the second équation of (4.28) we find 


As (x) r" + b 4 ( x ) r\ + d 3 ( x ) r, = o, 


(4.31) 


with 


d 3 (x) = è (</>' - 0 ) k 2 + ^k' 2 + e ! 

B 4 (x ) = ^ (0' — 0) L 3 + ^0 (. K 2 + L 3 ) + 6 (4.32) 

A (x) = jk<t>L 3 . 

Now it is necessary to search a common factor to reduce the degree of the polynomials 
A 5 , B 4 and D 3 . In fact, we hâve 


D 3 (x) = 0 (x) Pl (x), B4 (x) = 0 (x) Pl (x) and A 5 (x) = 0 (x) Pl (x) 


with 


$ = è 
ft = è 


— 20') 7T! + Ct/d] 

— 20') 7T + Jl 2 + 2o/0 + CU#' — 07r'] . 

Hence, we hâve (4.26) with deg0f < 3, deg0| < 2 and deg0| < 1. 
2) deg77 2 = 1, put H' 2 (x) = t. In this case, we hâve 

tr 0 = (H 2 r 0 y - H 2 r 0 . 


(4.33) 


Eliminating r 0 from the System (4.28), we obtain 


then 


h 2 (x) r' =r 1 (x)r ! + s 2 (x) r;, 
tr 0 = {k ' 2 - Rp r, + (k 2 + k ' 3 - s 2 ) r; + x 3 r". 


(4.34) 

(4.35) 


Using again the second équation of the System (4.28) to get 


where 


r } 1 (x) r" + p 2 (x) r; + rf (x) r 3 = o, 

V 3 = (0' — V 7 ) [| ( K 2 - Rl) H 2 + 07Ti - UJ0i\ = (0' - 0) /3g 
= (0 7 - 0) [f (A' 2 + L' 3 - R 2 ) H 2 + 07T - (U#] = (0' - 0) 0| 
v 1 = 7 (0' - 0) H 2 L 3 = (0' - 0) 0! 


4.3. Classical d-orthogonal polynomials 









Chapter 4. Characterizations of classical d-OPS 


58 


with deg /3l < 4, deg < 3 and deg P\ < 2. 

3) deg H 2 = 2. In this case, we hâve by using ( |4.34 ) 

H' 2 T 0 = (K 2 T 1 + LsT'J - R 1 T 1 - S 2 1^. 


Multiplying the first équation of the System (4.28) by H 2 H 2 we obtain the second 
équation of (4.26) where 


Pl = J (0' - 0) (K 2 - RJ H 2 + uH' 2 R x + 7r, H 2 H 2 
Pl = (J {P’ - 0) (K 2 + L' 3 - S 2 ) H 2 + üjH 2 S 2 + 7 tH 2 H' 2 
P! =ou'h 2 l 3 , 


with deg Pl < 5, deg P\< 4 and deg P\ < 3. 

For the first équation of (|4.26[), eliminate u\ from frorn (|4.28 ) we get 


/ii (x) Tx = m 2 (x) r 0 + n 3 (x) r' 0 , 


(4.36) 


where 

/ii = 9iti — tt 6 i, M 2 = 7 r (0' — 0) — 9u >', N 3 = np — c 06 , (4.37) 


with deg/ii < 1, deg M 2 < 2 and deg N 3 < 3. 

We distinguish two possible cases according to the degree of the polynomial /i 1 (x). 

1) deg/ii = 0, i.e. /i 1 (x) = C te := C. As in the first part, we obtain after simplification 
the first équation of (4.26) with 


a\ = 7Ti (0' — 0) — u [6 + M ' 2 
OL 3 — 0 7Ti — <jj 6 \ -\- M 2 -f- iVg 

»4 = iV 3 , 


where deg «4 < 3, dega| < 2 and deg < 1. 


2) deg/ii = 1, i.e. /i) (x) = k in the same way we hâve the first équation of (4.26) with 


a\ = üj'fJL 1 + |7Ti/ii + 7T0i 

«3 = 07/ii + |7Ti/ii (Af 2 + + 07Ti - U70i) + 7T (u;0i - 07Ti) 

«1 = \}i\N 3 . 


where deg al < 4, dega| < 3 and deg < 2. 


Conversely, it is enough to prove that we can deduce the System (4.28) from (4.27). 


Indeed, on account of (4.29)-(4.30) and (4.36 )-(4.37), the systern (4.27) may be written 


0 [w'r 0 + ttWi + Trr;] = oj [(0' - 0) r 0 + + evp , 

9 [TTifi + üüT 0 + W F 0 ] = TT [0iri + (0' - 0) r 0 + 0r / o ]. 


(4.38) 
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After adding uxj)T' 0 and 9nT\ respectively, to both sides of the first and second equality 
of (4.38), we obtain 


cf)X = uY 

ex = 7 tY, 


with 

À = [wT 0 + cuTq + ntT i + vrr;] 
y =[(0 , -‘0)r o + 0r{ ) + 0 1 r 1 + 0 , r 1 ]. 

Hence eliminating X and Y, respectively, we deduce 

(u9 — 7T0) Y = ( u9 — 7T0) X = 0, 
whence the resuit since u>9 — Tufi ^ 0. ■ 


4.4 Application 

Furthermore, in ra. the authors showed that there are exactly 2 d sets of d-symmetric d- 
classical polynomials (see also ra>- The rnethod used to obtain the differential équation 
of classical 2-OPS in [Tü] , can be used again for the d-symmetric d-classical polynomials 
to give the explicit form of solutions, fndeed, we hâve the System 


7n+d+l bn+1 ~ 7n+l] = 7n+l K+d+1 
in which we distinguish d+1 cases in eaeh case we hâve 


Tn+d+l] J 

) solutions, hence we hâve 

k ' 


d 


£ 



= 2 d 


solutions in total. In other words, solutions are 


S 0 ) 7°+i = 7n+i (trivial solution) 

s i) 7°n+ao = ( n n|p! o a ° ) ldn+a * and ^ dn + a = ldn+a for « ^ «o and 1 < a < d 

S 2 ) 7l +ao = ( ! SSf) 7 dn+a 0 an <* 7°n +Ql = ) 7dn+ Ql and 

7 dn+a = 7 dn+a for a ^ « 0 , oq and 1 < a < d 

Sd) 7 Sn+Q = 7 dn+a for 1 < a < d. 


4.4. Application 








Chapter 



Riordan Arrays and d-Orthogonality 


5.1 Introduction 

5.2 Riordan arrays 


The Riordan group introduced by Shapiro et al [72] is a set of infinité lower-triangular 
integer matrices where each matrix is defined by a pairs of formai power sériés g(z ) = 
Y^=o9n zU and f( z ) — fn z11 with go ^ 0 and f\ ^ 0. An infinité lower triangular 

matrix D = [d n ,k] n k >0 is called a Riordan array if its i th colnnm generating fmiction is 
g (x) [f {x)] 1 for i > 0 (the first colnnm being indexed by 0). With little loss of generality 
we also assume do,o = go — 1- The matrix corresponding to the pairs g,f is denoted by 
(g, f). One example of a Riordan matrix is the Pascal matrix 


P 


_i_ z —\ 

l-zl-z) 


10 0 0 
110 0 
12 10 


for which we hâve 



The group law is then given by 


\x 


1-z \l-z 


ig,f)-(h,l ) = {g(ho f),lo f). 


(5.1) 


The identity of this law is I = (1, x) and the inverse of (g, f ) is (g, /) _1 = (l/ (g o f) , f) 
where / is the reversion or compositional inverse of /. The reversion of / is the power 
sériés / such that (/ o /) ( x ) = x. We shall sometimes write this as f = Rev/. 
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A Riordan array of the form (g(x),x), where g(x) is the generating funetion of the 
sequence a n , is called the sequence array of the sequence a n . Its general terni is a n _fc. 
Such arrays are also called Appell arrays as they form the éléments of the so-called Appell 
subgronp. 

For a Riordan matrix D = [d n ,fc] n fc>0 , Rogers [67] lias fonnd that every element d n +i,fe+i 
can be expressed as a linear combination of the éléments in the preceding row starting 
from the preceding column. Merlini et al [BDJ hâve fonnd that every element in the 0 
column can be expressed as a linear combination of ail the éléments of the preceding row 
also see ra- That is, there exist unique sequences A = (ao, ai,...) and Z = (zo,zi,...) 
with «o 0, zo ^ 0 such that 

OO 

(1) ^n+l,fc+l hy d n .k \ J ■ {k)Tl 0, 1 , •••), 

j =0 

oo 

(2) dn+i,o zjdn t ji (ri 0,1,...). 

3=0 

The coefficients ao, ai,... and zo, Zi,... appearing in (1) and (2) are called the A- 
sequence and Z-sequence of the Riordan matrix D = (g(z), f(z)), respectively. Letting 
A (z) and Z (z) be the generating functions of the corresponding sequences, we hâve 

f(z) = zA(f(z)), 
g ( z ) = g (o) / (i — z (f (z))). 

We therefore deduce that 

A ( z ) = z /f( z ), 

and 

Z ( z ) = 1 /f( z ) (i -1/ (g ° f) (-)) • 

A conséquence of this is the following resuit wliich was originally established [5T] by Luzôn: 

Lemma 5.1 Let D = (g,f) be a Riordan array whose A-sequence, respectively Z- 
sequence hâve generating functions A(x) and Z(x). Then 

A — xZ x 
do,oA ’ A 

If (ao,ai,...) T is a column vector with generating function A (x), then multiplying 
D = (g, f) on the right by this column vector yields a column vector with generating 
function g(x)A(f(x)). 




5.2. Riordan arrays 
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5.3 Riordan arrays and <i-orthogonal polynomials 

Let {-Pn} n >o be a sequence of monic polynomials in V for which there exist complex 
sequences {/3 n } n>0 , (Xn,v), 0 < v < n such that 

P 0 (x) = 1, Pi (x) - x - An 

n— 1 (5.2) 

Pn+1 (x) = (X~ Ai+ 1 ) Pn (x) ~ X) Xn, V P V (x) , U > 0. 

The dual sequence {u n } n > 0 , «« e V of {P n } n > 0 is defined by (u n , P m ) := ô n>m , n,m> 0. 
In this case, we hâve 


Ai = (u n ,xP n (x)), n > 0, 

Xn,v = (u v ,xP n+ 1 (x)) , 0 < V < n. 


( 5 . 3 ) 


Let us consider d forms u 0 , ..., w^-i (d > 1). Let us now recall the following characteri- 
zation which we need in the sequel. 


Theorem 5.1 Let {P n } n>0 be a monic sequence of polynomials, then the following 
statements are équivalent. 

(a) The sequence {P n } n>0 is d-OPS with respect toU = (m 0 , •••, «d-i) T • 

(b) The sequence {Pn} n>0 satisfies a récurrence relation of order d + 1 (d > 1) : 

P-m+cL+l {x) (X fim+d) Pm+d {x) ^ ^ Xm+d-uPm+d—l—u {x) , TTl > 0 , ( 5 . 4 ) 

L — *is=0 

with the initial data 


P 0 ( x ) = 1 , Pi (x) = x - Ab 

P?n (x) = (x - An-l) Pm- 1 (x) ~ YZ =0 Tm-l-wPm-2-v (x) , 2 < Ul < d, 


( 5 . 5 ) 


and the regularity conditions 7^ + i A 0, m > 0. 


It is well known that we can express ( 5.4 ) as xP = JdP where P = (P 0 ( x ), Pi (x), 
and Jd = (ajj)°° =0 is a (d + 2)-banded lower Hessenberg semi-inhnite matrix, i.e., 


aij = 0, for j > i + 1 and i > j + d, 

1 A 2 > 0. 


( 5 . 6 ) 
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More precisely, the matrix Jd may be written 


/A) 1 

A 


J d = 


Y' 1 


Y ~ 2 Y - 1 


o 

î 

A. 


\ 


7i° 


il ■■■ 

7° 7s 


0 

1 


ld~ l 


Pd 


ld+1 Ad+1 


0 

1 


V 


/ 


and is called the monic Jacobi matrix of the monic d-OPS {P n } 


n> o- 


Now, when {P n } n > o i s n °t monic, then the linear récurrence (5.4) become 

yd 


Pn (x) = (a n X - fin- 1) Pn- 1 ( X ) - Yy Y-v P n-l-u (x) , U > 0. 

L - *V=l 


Write 


E TL 

a n,jX J 


then 


0 ^n^"ri— 1 , 0 ? 


^n^n— 1,0 fin—l&n— l,z—1 ^ ^ 7 /—1 r )n—v® i n—l i i—l—V’) 1 ^ ^ ^ 


n. 


(5.7) 


(5.8) 


(5.9) 


On the other hand, for a family of monic d-OPS given by (5.4), and by using (5.8), we 
can hâve 


n+1 n n —1 n—d 

^ ^ (•£ /3n) ^ ^ Tn ^ ^ ^n— l,k% ••• ^n—d+l ^ ^ Q j n—d,k% 

k =0 /c=0 fc=0 /c=0 

from which we deduce 


d -1 

®n+l,0 /dn®n,0 ^ ^ In—v ®n— 1— 11 ,0 (5.10) 

k=0 


and 


d—1 


®n+l,fc ®n,fc—1 f^n^n,k ^ ^ l n —v ®n—1 


— 1 —hh/c* 


k=0 


(5.11) 


We note that if /3 n and 7([ -î , 1 < i < d are constant, equal to f3 and 7 d_î , 1 < i < d, re- 
spectively, then the sequence (1, —/3, — 7 d_1 , —7 d_2 ,..., —7°, 0, 0,...) forrns an A-sequence 
of the coefficient array. The question immediately arises as to the conditions under which 
a Riordan array ( g , /) can be the coefficient array of a family of d-OPS. A partial answer 
is given by the following proposition. 


5.3. Riordan arrays and d-orthogonal polynomials 
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Proposition 5.1 [£$ Every Riordan array of the form 


x 


i+eïs«^'i+s:;v 


d -\-1 


(5.12) 


is the coefficient array of a family of monic d-OPS. 

Proof. The array [l/ (l + Y%=\ OkX k ^j , x/ (l + Ylk=\ OkX k ^ |4ÏÏJ lias a C'-sequence 

c{x) = Jf n >o c P x11 § iven b y 


X X 

1 + ES^" ~~ l-xC(x)' 


and thus 


d+1 

C{x) = ~'y' j 6 k x k ~ 1 . 

k= 1 


This means that the Riordan array ^1/ ^1 + Ylî =î OkX k ^J ,x/ ^1 + Y2Î=\ 
mined by the fact that 


is deter- 


®n+l,fc ®n,fc—1 T ^ ^ i,k ^ ■ h 0,1,2,... 

i>0 

where a n -1 = 0. In the case of ^1/ ^1 + Ylt=i @kX k ^j ,x/ ^1 + Ylt =î j j we hâve 

d 

®n+l,fc ®n,fc— 1 ^ ^ (5.13) 

v=0 

Working backwards, this now ensures that 

d 

En+lix') (*^ @l)En(,x') ^ ^ Qy+lPn—vi 

V=1 

where P n (x ) = J]Lo a re,fcZ n - ■ 

We note that in this case the (d+2)-term récurrence coefficients (3 n and 7^ _î , 1 < i < d, 
are constants. 

As an example, we hâve the following resuit 


Proposition 5.2 ffff The Riordan array (5.12) is the coefficient array of the modified 
d-orthogonal Chehyshev polynomials of the second kind given by 


P n (x) = U n (x- 0i), n = 0,1, 2,... 


5.3. Riordan arrays and d-orthogonal polynomials 
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Proof. We hâve bu ( 2 . 4 )] 


Thus 


Now 


--= y" u n (x)t n . 

1 oo 

Y J Un{x-e l )f 


1 — {x — 9f) t + Efc=l 9k+lt k+1 


n =0 


1 


i-(x-9 1 )t + z k =i d k+it k+1 V 1 + ELoW fc+1 i + ELcA+i tk+1 J 1 ~ xt 

Thus 


t 


1 + ELo °k+it k+1 i + ELo ^+i^ + 1 y l ~ xt 


Y J Un(x-9 l )V‘ 


72—0 


Also we hâve the following motivating resuit 


Proposition 5.3 fJfff Every Riordan array of the form 


1 ~ Efc=i ^kX k x 

i + Eiï^’i + Elïï^ 


is the coefficient array of a family of monic d-OPS. 


Proof. We hâve 


d -\-1 

B = (1 — A kX k , x) ■ A 

k =1 


where 


5 

A 



i-E 


d+1 

fc=l 


i+E 


d+1 
k =1 


1 


i+E 


d+1 

fc=l 


\kX k 
6f i x k ’ 


5 


_s_^ 

i+Efeiî y 

i+E2iî«fc*0 


7 


and (1 — Efii A kX k ,x ) is the array with éléments 


/ 1 0 0 0 

-À! 1 0 0 

—A 2 —Ai 1 0 

: —A 2 —Ai 1 


0 0 ... \ 

0 0 ... 

0 0 ... 


— Ki +1 

0 


—Ad + i 
0 


—A2 —Ai 


0 

1 


(5.14) 


(5.15) 


5.3. Riordan arrays and d-orthogonal polynomials 
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Now, by using the récurrence (5.13), we deduce 

bn+l,k ®n+l,fc n ^v+X^n—v,ki 

1 n @v+lQ"ri—v,k 


-Ai 


®n— l,fc—1 E.,-n ^v+l^n—1— 


-u./c 


- A 


< 2+1 


®n— 1— d,k— 1 E^O ^n+l®n— d— 1—■ 


î>,fc 


bn,k —1 E..-Q A^+l&n— v,k 


Hence, since b Uj k is expressed in terms of a ni fc, then as in proof of the previous propo¬ 
sition, il is a Riordan array of a farnily of d-OPS {Q n } , i.e. 


Qo{x) — 1 ) Qi(x) — x — 6\ — Ai, Qz^x) — x 2 — ( 20 i + Ai)æ + Ai$2 — A2 + 0 \ — 62, ■ ■ ■ 


Proposition 5.4 The Riordan array (5.14) ts the coefficient array of the modified d- 
orthogonal Faber polynomials given by 

P n (x) = F n (x - 9f ), n = 0,1, 2,... 

Proof. The d-orthogonal Faber polynomials {F n (£)} n>0 are studied by Douak and Ma¬ 
roni [32], and in another way by Ben Romdhane na. This family of polynomials is dehned 
with the help of the generating function as follows 


1 - Eh: v 


Tlius 


i.e., 


1 -xt + Eti 0k+it k+1 

1 - ES A k t k 

l-(x-9i)t + Eti 0k+it k+1 




n =0 


J^Fnix-ejr, 


n =0 


1 - ES +4 _ 

1 + Eh» 4+i ' 1 + Et» 4+i< l+1 ) 1 


t 


xt 


J2Fn(x-0 1 )t n . 


n =0 


Also note that, we can express the d-Faber polynomials in terms of co-recursive d-OPS. 
In fact, Ben Cheikh and Ben Romdhane m hâve shown that the cr-Appcll d-OPS are 
expressed in terms of d-orthogonal Chebyshev polynomials of second kind as follows 


<2+1 


Fn i,x) U n (x) T [ A kUn—k (.x) 


k =1 


5.3. Riordan arrays and d-orthogonal polynomials 
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In this case, since U Y (x) = U n (x), it follows that the Riordan array (5.14) is the 
coefficient array of the co-recursive d-orthogonal Chebyshev polynomials of second kind 
given by [69] 


d+1 


P n 


X 


= U n 


uc 


9\) + XkU n -k (x — Qi) 


k =1 


5.4 Stieltjes matrix 


Let L = (ln,k) nk>0 be Riordan and L be the matrix obtained from L by deleting the ûrst 
row, i.e., l nk = l n +i,k ■ For example, if / is the identity, then 

0 1 0 
0 0 1 

oo o i 


I = 


Observe that L = IL. There exists a unique matrix Sl = {snk ) nk >o such that LSl = L. 
That is, 

ln,k = y^ J Sjkl n -i,i, for n > 1. (5.16) 


i> 0 


We call this matrix the Stieltjes matrix of L. For example, the Stieltjes matrix corre- 


sponding to the Pascal matrix (5.1) is 


Sr = 


1 

1 

0 

0 

1 1 


0 

0 1 

1 

i.e., S L = 

= Sk L 


Theorem 5.2 $$ If L = ( g , /) is Riordan and Sl is (d + 2 )-banded lower matrix as in 


(5.1), then 


(a) Sl is in the foliowing form 

( A) 


S L = 


0 0 0 0 


\ 


7i _1 

it 2 Y- 1 


P i 

P 


7 d-2 y-1 


7i° : 

0 y 0 

: o 


ooo 
î o o 

P 


7 d-2 7 d-l ^ X 


(5.17) 


5.4. Stieltjes matrix 
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(b) f = x(l + f3f + j d 1 p + ...+'j°f d+1 ) 

(1 - xp 0 - 7 ? _ 1 x/ - ... - 7 %xf d y 


and 




Proof. Write Sl in the form (5.7) with C t (x) the generating function of the i th column 
of L, i > 0 . We hâve: Ci = gf l . By looking at the first column of LSl and L, we obtain 


C 0 ( X ) — 1 + ffiX + ff 2% 2 + ••• — 1 + ho^ "h 1 2 o ^ 2 + ... 


and 


Po — ho — Qi 
hoPo + 7i 1 = ho = 92 
hoPo + hi'yf 1 + 7i 2 = ho = 93 


hoPo + h\l\ 1 + i7i + 7? — ld+i,o — 9d+x- 

Multiplying the first équation by x, the second one by x 2 ,... and the last by x d+l and 
summing we get 

xp 0 g + X 7 d ~ l gf + ... + x^gf d = g- 1 . 

Whence g (x) = (l — x/3 0 — 7 f _1 x/ — ... — 7 ° xf d ) 1 . 

Repeat the procedure for z > 1, and in each case we try to find a System that contain 
ail coefficient P, ..., 7 °. One can see that we hâve 


sf = xg [r 1 + P,f + +.- + 7?+i f J+ '] , 

that is 

Ci = x (Ci -1 + PiCi + 7 î d + 1 1 C ',: + 2 + ... + ^ +l Ci + d) . 


Hence after simplification we find the expression of /. 

For instance, rewrite the equality (5.18) in ternis of j and subtract, we get 


d -1 


(Pi - Pi) / + £ ( 7 ; 


d—l—v 
i +1 


77;-”) / 2+ ”=0 


t?=0 


which is équivalent to 


(5.18) 

(5.19) 


Pi = Pi=P, lt'- v =l d - l - v , Q<v<d-\. 


This leads to the important corollary which is an analogue of 

Corollary 5.1 If L = ( g(x ), /(x)) is a Riordan array and P = Sl is (d + 2 )-banded 
lower given by (5.11), then L 1 is the coefficient array of the family of d-OPS. 


5.4. Stieltjes matrix 
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Proof. By Favard’s theorem, it sufïices to show that L 1 defines a family of polynomials 
{Pn( a;)} that obey a (d + 2)-terrn récurrence. Now L is lower-triangular and so L 1 is the 
coefficient array of a family of d-OPS P n (x), where 


( 1 ^ 


( Po(x) \ 

X 


Pl(x) 

x 2 

= 

P 2 (x) 

X 3 


P 3 (x) 

\ ; ) 


V : J 


We hâve 

S L ■ L- 1 = L- 1 L- L- 1 = L- 1 ■ 7 ■ L- L" 1 = L" 1 ■ ï. 

Thus 


Sl ■ L -(l,x,x,...) — L ■ I • ( 1 , x, x ,...) —L ■ (x,x , x 3 ,.. .) 


We therefore obtain 


( P o 1 0 0 0 

Y (- 1 p îoo 

Y ~ 2 Y -1 P i o 

; ryd-2 ryd—1 ‘ . 

7? : 7 d_1 P 

0 y 0 ••• 

v : o 


( Po(x ) 

\ 


( xP 0 (x) 

\ 

Pi(x) 



xPi(x) 


P 2 (x) 


= 

xP 2 (x) 


P 3 (x) 



xP 3 (x) 


V : 

) 


V i 

/ 


from which we infer that 


Pi O) =x-Pq, 

P 2 (x) = {x - P) Pi (x) - 7 i _1 


and 


p d+ 1 (x) =(x-P) P d {x) - Y~ l Pd -1 - -YPi Y) - 7i°. 

d—1 

P n +l( x ) = (x- P)Pn{x) -^Y^^Pn+d-l-vY), H > 1 . 

p=0 


Combining these resuit with previous proposition, we hâve 

Theorem 5.3 f A Riordan array L = (g(x), f(x)) is the inverse of the coefficient 
array of a family of d-OPS if and only if its Stieltjes matrix Sl is (d + 2 )-banded lower 
matrix. 


5.4. Stieltjes matrix 
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We give the following example 

Proposition 5.5 t 44\J The Stieltjes matrix of the inverse Riordan array 

1 x 

.i + Eti^’i + Et!^* 

left-multiplied by the k th binomial array 


x 


is given by 


1 — kx’ 1 — kx 


( @i + k 1 

@2 9\ + k 

9 3 


X 


1 — x' 1 — X 


P = 


0 

1 


0 \ 
0 


O 2 9i + k 1 
O3 d 2 


Od+l : 

0 9 d+ 1 

: 0 


Proof. We hâve 


x 


x 


1 + Efciî °kX k ’ 1 + E it\ 0kX k J V 1 + kx ’ 1 + kx 

1 


X 


1 + {9\ + k) x + Efc=i 9k+ix k+1 1 + {9\ + k) x + E&=i 9k+vx k+l 


And more generally, 

1-ESV 


X 


X 


1 + ES fet* ’ 1 + Et! <'*:V / V 1 + kx' 1 + kx 

1 - Et! Kx k 


X 


^ 1 + (9i + k) x + Efc=i 1 + ($1 + x + Efc=i @k+ix k+1 _ 

Therefore, the inverse of the last matrix has the Stieltjes array 


( 9i + k — Ai 
$2 — ^2 
93 — A 3 

9 d +1 — Ad+i 
0 


1 

9\ + k 


0 

1 


0 \ 
0 


02 9\ + k 1 

$3 9 2 

9d +1 
0 


(5.20) 


5.4. Stieltjes matrix 



























Chapter 5. Riordan Arrays and d-Orthogonality 


71 


5.5 Exponential Riordan arrays 


The exponential Riordan group is a set of infinité lower-triangnlar integer matrices, where 

2 

each matrix is defined by a pair of generating fnnctions g{x) — go + gi x + g 2 fr + ' ' ' and 
f(x) = fi x -\- fi\ ]- + ••• where go fi 0 and / 1 fi 0 . In what follows, we shall assume 

9o = fi = 1- 

The associated matrix is the matrix whose ?'-th column lias exponential generating function 
g(x)f {x) % /i\ (the first column being indexed by 0). The matrix corresponding to the pair 
/, g is denoted by [g, f], The group law is given by 

\9,f] ■ [M = [9(hof)Jof}. 

The identity for this law is I = [l,#] and the inverse of [<?,/] is [g, /] _1 = [1 /(g o /),/] 
where / is the compositional inverse of /. We use the notation elZ to dénoté this group. 

If M is the matrix [g,f], and u = (u n ) n >0 is an integer sequence with exponential 
generating function U (x), then the sequence Mu lias exponential generating function 
g(x)U(f(x)). Tlius the row sums of the array [g, f] hâve exponential generating function 
given by g(x)eA x l since the sequence 1,1,1,... has exponential generating function e x . 

We will use the following [25U2B], important resuit concerning matrices that are pro¬ 
duction matrices for exponential Riordan arrays. 


Proposition 5.6 Let A = (a n ,k) nk>0 = [ç{ x )i f{ x )] be an exponential Riordan array and 
let 



c{y) = 

C 0 + Cil/ + c 2 y 

2 +..., 

r(y) = r o + ny + r 2 y 2 + ... 

(5.21) 

be two formai power 

sériés that 







r(f( x )) = 

f( x ) 

(5.22) 




c(f( x )) = 

g'i x ) 
g( x ) ’ 

(5.23) 

Then 







(0 

(-1,0 

i 


(5.24) 


(«) 

^n+l : k 

r O a n,k-l + 

(ci-k + hri-k+Janj 

(5.25) 


i>k 


or, assuming Ck = 0 for k < 0 and rk = 0 for k < 0, 


^n+l,fc ^ ^ “1“ kri— 


i>k— 1 


(5.26) 


Conversely, starting from the sequences defined by (5.21), the infinité array (o, n ,k ) nk>0 
defined by (5.26) is an exponential Riordan array. 


5.5. Exponential Riordan arrays 
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A conséquence of this proposition is that the production matrix P = ( Pi,j)ij >0 for an 
exponential Riordan array obtained as in the proposition satisfies [25) [26] 

Pi,j = xÿ ( c i~j + jn-j+i) (c_ i = o). 

Furthermore, the bivariate exponential function 

(,b P (t,z ) = y^p n ,kt k ~i 
' n\ 

n,k 

of the matrix P is given by 


z) = e (c(z) + tr(z)). 


Note in particular that we hâve 


and 


r(x) = 


C[X = 


(5.27) 


(5.28) 


9(f(x)) 

Now, an important property of the matrix L which lias led to the proof of various 
results in this context is the following 

Proposition 5.7 ,11 1/ We hâve L = —L. 

dx 

Proof. Since 


d 

dx 


E: 

\n =0 


X 


n\ 


E 

\n =0 


X 


9n +1 ' 


n\ 


then equating the first columns of the matrices L and LS, we obtain 


Ai = ho 
hoPo + 7i 1 = ho 
hoflo + hilt 1 + 7i — ^30 


(5.29) 


IdoPo + Al7l 1 + ■■dd,d-ll\ + 7l —ld+l,0i 
and for n > d + 1, we see that one can lias 

InoPo + Al7l 1 + ■■■ln,d- l7l + 7lA,d = ln+ 1,0- 


(5.30) 


Multiplying the second équation of (5.29) by x, the third one by x 2 ,... and the last by x d 
and ( |5.30 ) by x n , and summing we get 

Ai (l + ^10^ + hox 2 + ...) + 7 i 1 { x + hix 2 + ...) + ... + 7° ( x d + ld+i,dX d+1 + ...) = g' (x), 


5.5. Exponential Riordan arrays 
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î.e. 

Po9 (®) + li 1 9 (x) f (x) + ... + 7 19 {x) [f (z)] d = g' 0) • 

Now, by equating the second columns of the same matrices we get 

1=1 
ho + Pi — hi 
ho + hiPi + 72 1 = hi 
ho + hiPi + ^3272 1 + 72 — hi 

ld +10 + ld+l,lPl + Wl,272 1 + ... + ld+ l,dl2 +72 = ^d+2,1 1 
and for n > d + 2, we see that one can lias 

ho + h,lPl + h,2l2 1 + ... + h,dl2 + h,d+ 17° = Wl,l) 


(5.31) 


(5.32) 


Multiplying the second équation of (5.31) by x, the third one by x 2 ,... and the last by x a 


and (5.32) by x n , and summing we get 

(1 + liox + I 2 ox 2 + ...) + Pi (x + hix 2 + ...) 

+72 _1 ( x 2 + l 3lX 3 + ...) + ...+ 72 ( X d+1 + l d+ 2,d+lX d+2 + ...) 

= {x + l 2 1 X 2 + ...) , 

i.e., 

g (x) + Pig (x) / (x) + 72 _1 g (x) [/ (x )] 2 + ... + 7 %g (x) [/ (x )] d+1 = g (x) /' (x) , 
whence 


1 + Pif (x) +7^ 1 [/ 0 x)Y + ... + 72 [f (x)] fl+i = /' (x) 


d +1 w 


Theorem 5.4 f If L — [g(x), /(x)] is an exponential Riordan array and Sl is ( d + 2)- 
banded lower matrix, then necessarily 


Sr = 


«o 

rir 1 


ho 


1 

01 1 


ho 2 hl 1 


hî 

h? 


0 

1 

«2 


0 

1 


hj- 1 


h2 


«d 1 

. ,d— 1 _ 

hd a d+l 


0 

1 


(5.33) 


■ 
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f d—i 

where {«&}*,>o is an arithmetic sequence with common différence a, and < 

1 < i < d, is an arithmetic sequence with common différence and 


k> 1 


ln{g)= I (a 0 + ^o 1 f(x) + ...+^ 0 [f{x)} a )dx, g( 0) = 1, (5.34) 


where ( a) k is the Pochhammer symbol and f is given by 


f' = i + af + ï> d ~7 <+1 , /(o) = o, 


(5.35) 


i=l 


and vice-versa. 

Proof. We consider the lower triangular matrix L with —g (x) [f {x)] k for the exponential 

rC! 

generating functions of the k th column, k > 0. We note that L is a Riordan matrix with 
exponential generating functions. 


From (5.35), we see that 


M g’ = Ai [x n ] g + Y. 7 "‘ K) gf\ 


2=1 


î.e. 


1,0 — fdo^nfi + 7l 1 ^n, 1 + ••• + 7l ln,d- 

Hence, for k > 1. we hâve 

1 \ ' i i 

— qf k ) =—n'f k + - fk-1 fl 

k\ yJ J k\ yJ (k- l)\ yJ J 

1 ( " " A d ^ " rh 1 1 + Æ/ + E 72 ~ 7 i+1 ) gf k ~\ 


= ^[Pof + E iVP ) g F + 


2=1 


(fc-1)! 


2=1 


then 


[F] ( ygf*) = [F] (u 1 ^, gf k 1 + [ x n } — \^ l \ f k 


k\ 


\x 


= \x 


(ij-' + klt 1 ) gfM + + (l1 tjï+pH gfM 


k\ 


k\ 


wF 9l " l+[xn] > k 

.d -1 .,0 


where 


whence 


-l [o? 71 ! — - - a f kJr ^ -i- \x n ~\ ————o 

+ [X 1 (k + l)\ 9J + - + [X l (k + d)\ 9J ’ 


Oik = A) + k/3 1 and /if * = (& + 1 )^ ( 7 ^ * + *) for 1 < i < d, 


ln+l,k ln,k— 1 C%kln,k Rfc ln, 1 "F “h 


(5.36) 


5.5. Exponential Riordan arrays 
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Finally, we easy observe that 


d.k — oik -1 = (3i : = a and 


,.d-i nd-i 

^ k ^ k ~ v = ~it l for 1 < i < d. 


(k + 1 ), (k) 


We give now an example for a d-OPS defined by a generating f on ction which can be 
also generated by a Riordan array. 

Let us consider a family of Hermite polynomials generated by the following generating 
function 


G (. x , t ) = exp < (x — b) t, — t 


* _ t i +1 


2=1 


i + 1 


£ 


n >0 Tl\ 


T Pn (X) t 


(5.37) 


This generating function satisfies the following first order differential équation 


G\ (x, t) = I x — b — dit 1 J G ( x , t ) 


(5.38) 


2=1 


Now replacing G ( x , t ) and G' ( x , t) by the right hand side of (5.37), we can easily obtain 


the linear récurrence satisûed by the sequence of polynomials {P n } n>0 as follows 

d 

Pn+ d +i(x)^(x-b)P n+d (x) — ^2a i {n + d) i P n+d - i (x), n > 0 (5.39) 


2=1 


which means that {P n } n > 0 is 2-OPS where (n + d) i — (n + d — i + l)j. 

By translation, we can take a = 0. In this case, we hâve G' x (x, t ) = tG ( x , t ), i.e., 


P'n+i (®) = (n + 1) Pn (x), 

that is, {P n } n> o is Appcll sequence, and then it is classical d-OPS. 
Now by choosing 

b = d and ai — i + l for 1 < i < d, 


(5.40) 


the corresponding Stieltjes matrix (5.7) may be written 
/ b 1 


ai 

b 

2a 2 

2a± 

d\üd 

6 a 2 

0 

ad (2 + d) 


0 


0 

1 

b 


\ 


0 0 

1 0 

da\ b 1 

d (d “h 1) o>2 (d d - 1) ®i b 


0 

1 
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where the common différence is 0 in ail case, i.e., 


tir 


t4-i 


(k + d) i (k + d-l) i 


= ai — ai = 0 . 


Therefore, from (5.34) and (5.35) we obtain 


/ (x) = x and ln(c/) = (b + Gpx + 2 a 2 x 2 ... + d\adX d )dx , 


i.e. 


g (x) = exp < bx + <5;X 


i +1 


2—1 


Corollary 5.2 [ffflfL = [g(x ), f(x )] is an exponential Riordan array and Sl in the form 
then Lis the coefficient array of the family of monic d-OPS. 

Clearly, L^ 1 is then the coefficient array of the family of polynomials {P n (x)}. Gath- 
ering these resnlts, we hâve 



Theorem 5.5 IfflJ An exponential Riordan array L = [g(x), f(x)] is the inverse of the 
coefficient array of a family of d-orthogonal polynomials if and only ifits production matrix 
P = Sl is in the form (5.33). 


Proposition 5.8 ffjfj Let L = [g(x), f(x )] be an exponential Riordan array with (d + 2)- 
banded lower Stieltjes matrix Sl- Then 

n\[x n ]g(x) (l, / (x),..., [/ (X)]*- 1 ) = U[x n ) := fi n , 

where U = (u 0 ,..., Ud-\) T is the linear functional that defines the corresponding family of 
d-OPS. 


Proof. Let L = (h,j)i,j> o- We hâve 


x n = Yj n ,iPi(x). 

i=0 

Applying 13, we get for each 0 < a < d — 1, 


u a (x n ) 


= u a (SiLo ln,iPi{x)) = Sr=0 l n,iU a {Pi(x)) 
= Y= 0 l n,A,a = L,a = n\[x n ]g(x) [f (x)]“ . 


Corollary 5.3 J Jf(Let L = [g(x),/(x)] be an exponential Riordan array with (d + 2)- 
banded lower Stieltjes matrix Sl- Then the moments /i n of the associated family of d-OPS 
are given by the terms of the first column of L. 


5.5. Exponential Riordan arrays 
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5.6 Sheffer Riordan array 


Let g (x) and f (x) defined as above. Then the polynomials P n (x) , n — 0, 1, ... defined 
by the generating function 

n>0 

are called Sheffer type polynomials with P 0 ( x ) = 1. 

Therefore, the set of ail Sheffer type polynomials {P n (x) = [x n ] g (x) e x ^ x ^} with an 
operation of nrnbral composition, forms a group called the Sheffer gronp. 

In the following paragraph, we give an example of Sheffer type Riordan array. Let 
us now consider a family of d-OPS of Charlier type defined by the following generating 
function ra 


G(x,t) = (l + ri) (l/T, exp 


i=l 


X > = 


V -,Pn(x)t r ‘ 

^n >0 n \ v ’ 


(5.41) 


This generating function is of Sheffer A-type zéro. Note also that when r —» 0, we get 
the previous example. 

The family of polynomials generated by (5.41) satisfy the following récurrence relation 

d 

P n +1 Qr) = (x + b - rn) P n (x) - y~] aj (n) i P n -j (x ), n>0 (5.42) 

i =1 


which rneans that {P n } n>0 is d-OPS where = (n + 1) 9 i+1 + n9i and (n) i 
In this case, the corresponding Stieltjes matrix take the following form 


( b 

1 

0 


0 

\ 

ai 

b — t 

1 

0 


0 

2a 2 

2 ai 

b — 2 t 

1 

0 


6 a 3 

6 a 2 

3a\ 

b — 3r 




(n - i + 1).. 


V 


d\dd 

0 dd( 2 + d) d 

0 


da i b — dr 

d (d T 1) 0/2 (d T 1) cp 


/ 


r d ~ 

where the common différence of {o^ = b} k>0 is r and it is 0 for < 
Therefore, from (5.34) and (5.35) we obtain 


i ) k >0 


f(x) 

Mg) 


= b + xe TX and 


— J b + ai (b + xe Tæ ) + 2a 2 (b + xe TX ) 2 + ... + d!a d (6 + xe TX ) c 


dx, 


5.6. Sheffer Riordan array 
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i.e., 

d 

ln g ( x) = E ô t x l exp {irx} . 
2=1 


5.6. Sheffer Riordan array 



îhapter 


6 


Conclusions and future directions 


In this thesis we are concerned with Riordan arrays and ci-orthogonal polynomials. After 
giving sonie preliminaries results in the two first Chapters we focused in Chapter 3, 
on the construction of the inner products which generate the sequence of d-orthogonal 
polynomials in the sense of Sobolev. Our perspective in the future is to study d-classical 
of Sobolev and d-serni classical orthogonal polynomials of Sobolev. 

In chapter 4, we characterize the d— classical orthogonal polynomials. . We présent an 
algebraic theory of classical d—orthogonal polynomials and we want to hll in sonie gaps. 
We broaden and close sonie inclusions that exist and are known perhaps as conséquences. 
In the future we aim to extend the combinatorial theory of orthogonal polynomials created 
by Vienno to d-orthogonal polynomials. Indeed, a second method of calculating moments 
lias been discovered by using a non-trivial combinatorial model from a generalization of 
Pollaczek polynomials that can be generalized to ail Appcll polynomials. 

In Chapter 5, we pioneer the study of d-orthogonal polynomials and Riordan arrays. 
More precisely, we hâve extended the results of Riordan arrays and orthogonal polynomials 
in such a way to give a signihcant generalization from the orthogonal case to the d- 
orthogonal case. We show that the Riordan array, in the constant coefficient case, is a 
good vehicle for such studies. A principal discovery is that the d-orthogonal polynomials 
are characterized by the Rordan arrays. Our perspectives are then to extend this idea to 
algebraic structures in d-Hankel matrices and d- Hankel-plus-Toeplitz matrices relating 
to eOclassical orthogonal polynomials. We aim also to extend our research to d —continued 
fraction and d-semi classical orthogonal polynomials. 
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